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COHEIGUPATION OPTIMIZATION OP TRMSMISSION LINE TOlOTS 

by 

Eaxuppaiab Eumarasaray 

nth the enormous increase in the amount of electrical 
energy transmitted and the high voltage of transmission, the 
shape and size of the transmission line towers have received 
considerable attention in recent years. Since the tower 
structure is repetitive in nature, any saving in weight will 
build up the saving of material cumulatively. Hence the 
tower is most suited as an optimum design problem, 

V 

Sixties witnessed enough success in optimum structural 
design. Till early seventies these studies were by and large 
restricted to a fixed geometry of the structure. Off late 
attention is rightly focussed towards configuration optimization. 
In fact evolution of an optimal configuration should be the 
first step in any attempt of structural optimization. This 
thesis is directed towards the study of optimum configuration 
of Transmission Line Towers. 

Automated optimum design of transmission line towers, 
modelled as a space truss, under static and d5mamic behaviour 
is Carried out with the purpose of generating optimum 
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conf iguration. Tlie objecttve in the present work Is to 
minimize the total weight of materials used including the 
weight of secondary members. The weight minimization is 
carried out subject to the limitations on static and dynamic 
stresses in the individual members and the req^uirement of 
compatibility of the tower configuration. The design 
variables chosen are the geometrical dimensions of the 
transmission line tower, viz,, base width of tower and the 
panel heights in the body of the tower. 

The justification to model the tower as a space 
truss both for static and dynamic analysis has been outlined 
in the present work. The decision to model the tower as 
space truss for dynamic analysis has been arrived from the 
study of eigenvalues of the tower. The configuration of 
the tower is generated automatically. Study of member forces 
under all possible independent loading conditions reveals 
that only three loading conditions are critical. The static 
analysis is carried out considering the maximum wind load as 
a static load. In dynamic analysis, wind load is considered 
time dependent. Dynamic analysis is also carried out for 
earthquake loading. 

The optimum design problem is formulated as an 
unconstrained minimization problem, jill the constraints 
are handled implicitly. This is necessitated since the 
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objective function is not an explicit function of design 
variables, viz,, base width and panel heights. Powell's 
method has been used which turns out to be the obvious choice 
for seeking the solution of such unconstrained minimization 
problem. 

A study has been conducted to assess the effectiveness 
of the two methods of analysis, viz. , stiffness method and 
flexibility method, to be incorporated in optimum design. 
Stiffness method has been found to be efficient from computer 
time and storage requirements. The design of members of the 
tower is carried out as per the code IS: 802-1973, using 
available angle sections. , Prom the areas picked up for the 
individual members, it is observed that certain angles are 
always left out during selection. This led to the conclusion 
that some of the available angles are not rational. Prom 
this observation a study has been conducted to bring forth 
a set of rational angles. These are presented in i^pendix 2 
of the thesis. Using these rational angle sections, the 
transmission line towers are again optimized for static 
loading which gives further reduction in weight. 

Dynamic analysis of the transmission line tower is 
carried out through modal superposition. Dumped mass 
formulation is used in the present work. The resulting 
eigenvalue problem of the space tower is solved through 
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SubspaCG iteration metlaocL, The eigeasolution is a series 
of bonding modes and torsional modes in the lowest end of 
the spectrum* A parametric study has been carried out to 
study the variation of eigensolution with the height of 
the tower. Prom the results of parametric study, 
simplification on the modelling of slender towers for 
approximate dynamic analysis is recommended. The results 
of dynamic analysis, carried out separately for wind and 
earthquake loads, are tabulated. 

The optimization study is carried out for 15m, 20m 
and 25m body heights of transmission line towers under static 
and time dependent loading cases. The effect of damping on 
the optimum weight of tower is also studied. 

Some of the salient conclusions of the present 
study are: (1) the configuration corresponding to optimum 
tower is more sensitive to dynamic loads compared to the one 
obtained by considering only static behaviour; (2) for the 
design of transmission line towers, dynamic analysis under 
wind load should be carried out to have more realistic 
response prediction than considering equivalent statie wind 
loads; (3) the conventional model of plane frame for the 
dynamic analysis of slender towers is confirmed to be 
reasonably accurate through the present study; (4.) 'earthquake 
load is not critical for the design of transmission line 
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tourers; and (5) availalDle Indian Standard angle sections 
can "be revised to more rational sections from the view- 
point of minimum weight design. 
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HTTRODUCTIOIT 


1*1 Generajl 

The electric power generated at the themal, hy^iro 
and nuclear power plants is distributed far and wide by a 
network of transmission lines. Transmission lines can, 
therefore, be compared to the circulatory system in the 

I 

human body which distributes the energy required by the 
various parts of the body. Transmission lines are, as of 
today, a set of overhead conductors and a groundwire which 
transmit the electrical energy as high voltage current. 
Supporting structures are constructed at intervals to 
keep these lines at a clear height from the ground level. 
These structures are known as transmission poles or towers. 

Transmission poles are used for relatively low 

voltage transmission (generally uptp 152 kilovolts) and 

are generally made of wood| rolled or built up steel sections 

or prestressed concrete. Transmission towers are generally 

made of mild steel. Sometimes transmission towers of 
1 2 

aluminum and reinfoiced concrete are used for transmission 
of i^oltagpe abdve 152 KV, 

The transmission stfuotures can be divided into two 
Categories depending on their structural action, namely, 



2 


guyed or self-supporting, -n guyed towers, guy wires 
caxry most of the transverse loading and the central mast 
carries the vertical loading, A self supporting structure, 
as the name implies, carries the applied load all "by itself. 

.The guyed tower * requires less structural steel 
but more right of way and often more construction effort. 

The tower can be completely assembled on the ground and 
erected as a whole with the help of a crane. The self 
supporting tower is heavier, but has the advantage of 
occupying less ground area and providing greater stability, 

Guy anchors require excellent compaction to avoid yielding 
and slackening of the guys, whereas the average self 
supporting towers can tolerate a certain amount of substandard 
backfilling of its footings. Also, the guyed towers are 
vulnerable to sabotaging. Self supporting towers are, 
therefore, more widely used than the guyed towers. Hence 
the present study is directed towards the self supporting 
transmission towers of the kind shown in Fig. 1,1, 

Transmission towers are divided into three main 
Categories depending on their location. These are: 

(a) Tangent or suspension towers (Fig. 1.1): These are 
used in the straight sections of the lines. They 
are designed to support the wires when subjected 
to wind and ice loads and frequently are designed 




Typical double circuit transmission 
line tower 
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for a small angle cl ange in tlie line. They are also 
designed for the unbalanced pull due to a broken 
wire on one side of the tower, 

(b) Angle towers: These are located at points where 
the line changes direction. They are designed to 
carry the same load that are carried by the tangent 
tower in addition to a transverse pull due to the 
change in the line direction. In order to utilise 
the structure most efficiently, the structures are 
usually placed so that one of the axes of the cross 
section bisects the angle formed by the conductors 
as shown in Pig, 1,2, 

(c) Bead end or anchor towers: These are designed to take 
the dead end pulls from all the wires on one side, 
together with wind and vertical loads. 

The trend in the electrical transmission technology 

is towards the use of veiy high voltage for transmission. 

Voltage of the order of 750 are already popular in 

8 

the western countries. Billard and Hileman have reported 
that transmission of 1100 'PV is technically feasible. This 
increases the load on the structure and also the clearance 
requirements and consequently larger structures are required. 
An accurate analysis of the tower under these heavy loading 
and a synthesis of the analysis and design procedure to 
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Fig.1'2 Classification of towers based on location 


economise the cost of the "towers is warranted. The 
oomp"uterised analysis and optimum design is the natural 
choice of the means to achieve this end. 


6 


1,2 Review of Literature 

The first known paper on transmission line towers 

Q 

is hy Bergstron, Arena and Kramer , This is the state of 
art type paper, describing the salient fea’tures of the 
design of self supporting transmission towers subject to 
the then existing code of practice — the Fational Electric 
Safety Code, the shortcomings of this code and the necessity 
to modify the code, (The main points suggested in the 
paper are as follows; 

1, Selection of design loading conditions comparable 
with the anticipa,ted loads rather than provision 
for every remote contingency. 

2, Tillillingness to consider rectangular tower cross 
section with its increased complexities, 

3, Extensive studies on various configuration to achieve 
optimum strength with least material, 

4, Careful checking and rechecking of mraibers to ensure 
selection of minimum sizes and the use of a computer 
program to select members automatically. 

Incorporation of high strength steel in the design 
at advantageous points. 


5 
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The method of least work has 16011 used in this paper for 
the analysis of the tower as a space truss. 

Selection of site and heights of towers, depending 

upon the terrain, is a major factor in the design of 

transmission system. Dynamic programming or recursive 

optimization is especially useful for multistage process 

in which decisions taken at one stage do not affect the 

previous stages. Dynamic programming approach has been 

10-12 

successfully employed to design the transmission line 

tower location along the route so that the construction 
costs are minimized, while the physical constraints are 
satisfied. The input data required include the survey 
data of the route and the choice of available towers of 
suspension type and of angle towers. 

Having fixed the tower location along the route 

of the line, the next job is to decide the outline diagram, 

length of crossarm and conductor spacings, tower width at 

13 

the base and top hamper, and type of bracing pattern . 

Gopalan^^ has presented the procedure for the design of 

barrel towers and shown the transverse and longitudinal 

faces of 132 TY double- circuit barrel type tower. The 

13 

paper by Do, et al describes the ccmputer program 
developed for the analysis and design of steel transmission 
towers or general truss. The program capabilities consists 
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of (1) generating S57mmetric portions of geometry j 

(2) identifying and stabilizing planar (unbraced) joints; 

(3) handling only tension members; (4-) accounting for the 
buckling of members by using a bilinear force- displacement 
relationship; and (5) varying some hey joints to optimize 
the tower geometry. 

In the western countries, computers have been used 

to design completly an entire stretch of transmission line, 

which include the selection of conductor t 5 rpe, layout of 

the towers along the line, fixing the dimensions of the 

towers to satisfy clearence requirements, design of the 

towers and inventory control of structural and electrical 

components, in excell^t account of an integrated 

transmission line system design developed by Ms. Sargent 

1 6 

and Lundy Engineers is given by Beck . The fundamentals 

involved in the design of a transmission tower have been 

17 

described vividly by Marjerrison , 

With the introduction of extra high voltage (EHY) 

for transmission of electrical power, and the availability 

of the digital computers for analysis and design, a rational 

assessment of loads, allowable stresses and overload factors 

have become necessary and this was clearly brought out by 
is 

Parr, et al . 1 review of past practices in the calculation 

of loads on transmission line towers, as used in the design, 
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is given and a new concept in design is suggested using 

1 9 

more realistic ultimate load approach.* Campbell has 
presented a non-linear analysis of an elastic transmission 
line for conditions involving imbalance in conductor 
tensions. Stiffness properties of structures, conductors 
and insulators are derived considering structure and 
insulator displacements. 

Early design of transmission line towers were 

carried out by analysing the approximate determinate 

structure either by resolving the joint forces or by 

graphical methods. Today, with the advent of high speed 

digital coijiputers, it is possible to analyse the entire 

20 

tower as an indeterminate structure, Natarajan has 

presented a simpler version of the frontal solution 

technique suitable for structural analysis with large 

number of members. The method has been used to analyse a 

transmission tower and the resulting stress distribution 

in the main members are discussed. The elastic transmission 

system power flow equations bear a striking similarity to 

force- displacement characteristics of pin- joint structures 

such as electrical transmission towers. This similarity 

is successfully exploited in order to radically increase the 

21 

speed of anal3rsis of transmission line structures 
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In general transmisi-ion line towers are tested, 

tinder the expected loads, so as to guarantee the safety 

of the system under v^orking conditions, Murthy and 
2 ? 

Mukerjee have reviewed the methods and practices followed 
in some of the advanced countries as well as in India as 
regards to the testing of transmission line towers. A 
uniform procedure for testing the steel transmission towers 
is presented with suggestions for the improvement of the 
testing practices. 

23 

The paper hy Anaston describes a computer program 
written to design optimally a complete transmission tower 
of given basic dimensions. The tower, divided into five 
parts, has been analysed as a statically determinate system 
and the members are designed according to code specifications. 
The variables involved in the procedure are the width of 
the basket, width of the base, number of panels in the basket 
and number of panels in the body. The input consists of 
basic dimensions and loads. The configuration has been 
generated by the computer and the analysis and design are 
performed. The output consists of the values of the variables 
mentioned above, and the sizes and lengths of all the members, 
A distinguishing feature in this study is that the weight 
of the tower, excluding the weight of secondary members but 
including that of anchors, is minimized. The paper, however, 
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does not mention how the vaxdahles were varied and the 
amount of variation in each variable . 

Raj and hurrant^'^ have presented dynamic programming 
foimulation for structural design optimization in general 
and discussed possible extension to transmission line tower 
configuration, 

25 

Verma has analysed basket portion and body of 

the tower separately and developed a program, using force 

method of structural analysis, for minimum weight design 

of transmission line tower with square base. The design 

variables y/ere panel heights of the body and base width of 

the tower. Using Univariate method of unconstrained 

minimization, he concluded that panel heights have considerable 

effect on the minimum weight of transmission line towers, 

26 

Kuzamanovic, et al have reported considerable weight 
reduction by using three legs instead of classical square 
base for the transmission line towers, 

Uo effort seem to have been put, so far, towards the 
development of rational optimum design of transmission 
line towers subject to both static and time dependent 
loading, 

1,2,1 Optimization of Truss»-Iike Structures 
27 

Michell ' has solved the following problem of 
stiructural optimization: G-iven an equilibrium systan of 



12 


forces (loads and reactions), to design a plane truss with 
bars of a given e'^.astic material in which the given forces 
produce axial stresses within a given allowable range, 
with the amoimt of material required for the bars is as 
small as possible, Michell's optimal structures, however, 
are not trusses in the customary sense of the term but 
truss-like continua with an infinity of joints. Structural 
optimization problem today is nothing different from that 
of Michell’s , but stated in the generalised form. The 
prublem statement varies only depending upon the objective 
of the problem, 

28 

Excellent surveys by WasiutynsldL and Brandt in 

2Q 

1963 and Shcn and Prager ^ in 1968 provide a complete 

historical development of the analytical methods of 

structural optimization. Mathematical programming techniques 

applied to structural optimization problems were reviewed 

50 

by Kapoor and Hariharan till 1972, Recently Rozvany and 

31 32 

Mroz ^ have reviewed the progress in structural optimization 
at a conceptual level. 

Early works on structural synthesis were applied to 
truss-like structures, being simplest of stiuctural forms, 
to arrive at a minimum weight. The attention was focussed 
to get minimum cross section for the members of the truss 
subject to non-zero value of cross section, stress and 
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local instability constraints ♦ Th.is resulted in a set 

of necessary and sufficient conditions for local optimality 

54 

of a fully stressed design under a single system of loads 
35 

lev has presented ah algorithm that starts with an 

indeterminate truss and at each cycle of the procedure 

one bar (or more) is eliminated, thus yielding an optimal 

truss with a lower degree of indeterminacy, until the final 

solution is reached; Based on the least-weight indeterminate 

and least- weight determinate tiruss, useful bounds on the 

truss weight has been obtained which can be used as a 

36 

criteria for efficient design; Prager considered the 
near optimal trusses, The upper and lower bound optimal 
layout of truss was discussed that requires an amount of 
material exceeding the minimum possible amount by at most 
a known percentage. 

With the initial success in minimum weight design 

of simple structures, the question arose whether minimum 

weight design and fully stressed design are the same, 

37 

Dayaratnam and Patnaik have reported that the fully 
stressed design of a statically indeterminate truss is not 
feasible for single load condition, There exists a 
relationship between the number of load conditions and the 
order of indeterminacy of the structure for a feasible fully 
stressed design. 
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The general structurcJL optimization problem is to 

get tbe minimum of the objective function under a set of 

constraint functions which are, in general, nonlinear. 

These nonlinear functions of constraints and objective 

function can be expanded in Taylor’s series about the 

initial design point. If the series is approximated by 

retaining only the first two terms (constant and linear 

term) , then the problem becomes a linear programming problem 

which is comparatively easier to solve. The solution can 
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be obtained in the direct form or in the dual form 

depending on the computational effort req.uir8d. In general 

this t37pe of transformed SIP (sequence of linear program) 

problems produce unacceptable designs. To eliminate 

59 

infeasible designs, a set of pseudo constraints are added 
to ihe original linearized co-nstraints. This way the problem 
size becomes large which requires more computational effort 
and slow convergence to the optimum of the original problem. 

Based on the analysis of direct (mathematical 
programming) and optimality criteria (indirect search) 
methods for structural optimization, efficient hybrid methods 
were suggested exploiting common features of both methods, 
Pleury and G-eradin^ have shown, through numerical applications, 

significant improvement in convergence properties with hybrid 

41 

method. Arora and Haug have discussed the application of 
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hybrid methods to dynamic response problems, 

With the success in application of mathematical 

programming methods to statically loaded structural 

optimization problems, attempts were made to optimize 

structures under dynamic and stability constraints. Fox 
42 

and Kapoor have optimized the truss- frame structures, 
using the method of feasible directions, imder dynamic 
response regime. The members were designed under peak 
stresses. Cassis and Schmit have used the first-order 
Taylor series esipansion of the dynamic response functions 
with respect to the design variables. The Davidon-Fletcher- 
Powell algorithm has been used for the unconstrained 
minimization, 

Khot, et al^"^ considered the minimum weight design 

problem with stability constraints. The method was programmed 
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for trusses and frames, Schmit and Ramanathan have reported 
the minimum weight structural design of truss and wing 
structures in a multilevel approach. The overall proportioning 
of the stiucture was achieved at the system level subject to 
strength, displacement and system buckling constraints 
while the detailed component designs were carried out 
separately at the component level satisfying local buckling 
constraints, Kie systaa buckling behaviour was handled 
by incorporating the geometric stiffness matrix capability. 



Sctaiii; and have developed programmes 

that combines finite element analysis techniques and 
mathematical programming algorithm using an innovative 
collection of approximation concepts. Design variable 
linking, constraint deletion technique, and approximate 
analysis methods were used, Gkieir earlier paper^ discusses 
the application of AGGESS-1 ( J^proximation Goncept Gode 
for Efficient Structural Synthesis) to systems represented 
by truss, constant strain triangle and symmetric shear 
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panel finite elements. The later advanced version ACGBSS>-2^' 
includes problems involving fibre composite structure, 
thermal effects and natural frequency constraints in 
addition to the usual static stress and displacement 
limitations. 

Till early seventies, studies on structural 
optimization were, by and large, restricted to a fixed 
geometry of the structure, Of^ late attention is rightly 
focussed towards configuration optimization. In fact 
evoluation of an optimal configuration should be the 
first step in any attempt of structural optimization, 

48 

Tanderplaats and Moses^ have reported significant 
weight reduction when geometric changes are included in 
the optimization process. The member areas and joint 
coordinates were treated as design variables. The design 
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problem was, in effect, divided into two separate, but 

dependent, desigi' spaces one for member sizes and one 

for coordinates. Ihile changing coordinate variables the 
member areas were treated as dependent variables. Fully 
stressed design was assumed optimal. Eunar and Ohan^^ 
have formulated the minimiun weight design of structures 
with variable geometry. The member sizes were esipressed 
in terms of configurational variables satisfying the 
imposed constraints by use of optimality criteria. The 
problem was solved via an unconstrained minimization 
algorithm, 
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Bennett*^ has presented a two step procedure in 

which the first step include only the cross section variables 

and for the second step both cross section and layout 

variables were used. Method of feasible directions has 

been used for the minimization. Illustrative examples 

5 1 

include frames, lipson and igarwal have applied complex 
method of Box to optimize planar truss structures. The 
design variables included geometry and member area using 
the ■ available sections. This method has been extended 
successfully to space trusses and comparison was made 

with minimum weight design obtained with continuous variables. 
Favourable design improvements and rates of convergence 
were demonstrated. 
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P resent Work 

The computational capability for the automated 
optimum design of transmission line tower modelled as a 
space truss under linear elastic static, and d 3 mamic 
behaviour is developed in the present work with the purpose 
of generating optimum configuration. The objective function 
chosen is the total weight of materials used including the 
weight of secondary members. The weight minimization is 
carried out subject to the limitations on static , and 
dynamic stresses in the individual members and the requirement 
of compatibility of the tower configuration. The design 
variables chosen are the goemetrical dimensions of the 
transmission tower, viz, , base width of tower and the panel 
heights of body of the tower. 

The necessity of modeling the transmission tower 
'as a space truss structure both for static and dynamic 
analysis is brought out in Chapter '2. Automated generation 
of tower configuration is the first step towards the study 
of configuration optimization. This is achieved in two 
parts. The first part generates the fixed data such as the 
joint connectivity of members, and type of members and 
this is done only once during the whole optimization routine. 
This subroutine also generates the coordinates of the joints 
in the basket portion, as this part of the tower cbnfi^ration 
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is fixed. The second part generates the variable data 
such as the coordinates of joints in the body of the tower 
ahd length of secondary members, which get modified at 
each design point. This is described in CShapter 2. 
Transmission line tower is to be analysed for a set of 
load conditions as per IS:802-1967^^* These are discussed 
in Chapter 2, Out of all possible conditions only a few 
turn out to be critical load conditions from the view 
point of design. These are brought out in Chapter 2, 

In the 1973 revision of 13:802, the wind pressure values 
have been increased by 30 percent to account for d 3 mamic 
effect, irrespective of the height of tower. The constant 
increase in wind pressure for all heights does not seem to 
be correct, iny measure to account for the djmatnic effect 
has to bring- in the stiffness of the structure. . This is 
important while designing slender structures like transmiss-’ on 
line towers. Hence a deteministic d 3 aiamic load has been 
approximated wherein the maximum wind pressure at various 
heights correspond to those values given in 13:802-1967, 

This is discussed in Chapter 2, Earthquake loads for the 
design of transmission line tower are also discussed in 
this chapter. 

Chapter 5 deals with the formulation and method 
of Optimization of the resulting optimum configuration 
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design problem. The implicit nature of the objective 
function has made it necessary to formulate it as an 
unconstrained minimization problem, Powell's method 
turns out to be the obvious choice for seeking the solution 
of such an unconstrained minimization problem. 

Structural analysis consumes a major part of 

computer time in the design of transmission line tower. 

Optimization process being a repetitive analysis design 

cycle, it is of utmost importance to use an analysis 

technique which consumes less time. Hence for the present 

work both flexibility and stiffness methods of stiuctural 

analysis are studied to assess the effectiveness of one 

over the other from the computational point of view. 

Stiffness method has been found to be efficient from 

computer time and storage requirements. This is discussed 

in Chapter 4. The design of members of the transmission 
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tower is carried out as per 13:802-1973 using available 
56 

an^e sections . The area and other sectional properties 
of available sections are stored in the memory of the computer 
in ascending order of their axeas. The results of optimization 
of the tower configuration under static loading is presented 
in Chapter 4« Prom the areas picked up for the individual 
members, it is observed that certain angles are always 
left out during selection. This led to the conclusion that 
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some of the available angles are not rational. From this 
observation a study has been conducted to bring forth a 
set of rational angles. These are proposed and presented 
in ippendix 2, Using these rational angle sections, the 
transmission line towers are, again, optimized for static 
loading which gives further reduction in weight, OOhe 
details are presented in Chapter 4. 

Chapter 5 deals with dynamic analysis of the 
transmission line tower. Modal superposition method, 
used in the present work, is briefly discussed in this 
chapter. !IIhe resulting eigenvalue problem is solved 
using subspace iteration method. A parametric study on 
eigenvalues by varying the height of tower is carried out. 
The results of dynamic analysis for wind and earthquake 
loads axe presented separately. The optimum design under 
static loading condition, obtained in Chapter 4, is taken 
as starting design for carrying out optimization study 
under dynamic loads. These results are also presented 
in Chapter 5. 

The conclusions drawn from the present study 
and suggestions for further extension of the present work 
axe presented in Chapter 6. 
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TRMSIilSSION LINE TOlER STRUCTURE 

2,1 In troduction 

The transmission line tower consists of two parts 
(1) the basket or cage, and (2) body as shown in Eig. 1.1. 

The basket portion of the tower has the cross ams 
projecting from the tower shaft which in turn support the 
suspended conductors. The vertical spacing of these cross 
arms mainly depends upon insultation requirements.- This 
portion of the tower is detailed from electrical considerations. 

The body of the tower extends from the lowest cross 
arm upto the support , The height of body of the tower may 
vary from tower to tower depending upon the local geographical 
conditions. The configuration is decided by the width of 

base, number of panels and height of each panel , 

/ 

The tower structure consists of one dimensional 
slender members, as is seen from Eig* 1.1. These members 
can be classified as leg members, horizontals and bracings. 

They are generally ri vetted at ends to form joints . Sometimes 
the joints are also bolted connections. These type of 
connections lead to partially fixed condition at the joints. 
However, the slender nature of the members of the tower make 
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the effect of partial fixity of the joints negligible. 

Hence the joints of the tower are assumed to be hinged for 
the purposes of analysis. This idealization results in 
axial forces of tension or compression in the members 
when the tower is subjected to joint loads. 

Haring identified the structure, deciding a suitable 
model for analysing the tower structure, subject to static 
and dynamic loading, is the first step in the present 
study. This is discussed in Section 2,2, Automated 
generation of tower configuration as well as design load 
conditions are also discussed in the subsequent sections. 

2 *2 Mo deling of the T ower Structure 

In the design offices, the analysis of the tower 
structure is oversimplified. The simplification is two- 
fold: first, the tower which in fact is a space structure, 
is idealized to a planar truss and secondly the analysis 
of planar truss does not take into account the presence of 
redundants. In other words planar truss is analysed as a 
determinate structure. These simplifications are not 
justificable for the following reasons: 

(1) Projected dimensions of the tower, on a vertical 

plane, are considered. 
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(2) Loading on 1:110 tov;er, particularly the situations 

in which one of the conductors is broken, is not 
planar. 

Moreover, the simplifications discussed above which were 
primarily due to computational ease are no more warranted 
with the availability of high speed digital computers. 

A study has been conducted by Sondth and 
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Mukhopadhyay , modeling the tower as planar as well as 
space truss subjected to static loads. The computerised 
results obtained for the member forces show that, in 
general, the forces are overestimated in the planar model. 
Moreover, some members also seem to be highly understressed. 
Hence the tower structure has to be modelled as a space 
truss , 

2,2,1 Mod el for Static Analysis 

Modeling the transmission line tower as space 
truss result in three degrees of freedom at each joint, 
one each along the coordinate axes. The base support is 
assumed fixed for translations* Hence the number of degrees 
of freedom for a transmission tower with n joints (excluding 
supports) is 3n, 

2 *2 .2 Model for Dynamic Analysis 

The transmission line towers vibrate transversely 
in the bending modes, when all wires are intact. However, 
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when ground wire is broken the tower vibrates in bending 
in the longitudinal direction of the line. ]?urtheniiore, 
when one of the conductors is broken, the tower vibrates 
in torsional modes. These considerations make it necessary 
that the transmission tower should be modelled as a space 
frame . 

The main reason to model the stinicture as frame 

is that the individual member frequency of long member, 

assumed as simply supported at ends, may have fundamental 
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natural frequency less than that of the structure . In 
such a situation that particular member may vibrate with 
high amplitude and lead to failure which in effect mean 
collapse of the structure. 

The members of the transmission line tower are 
very slender. Hence the rotational resistance offered by 
the members at the joints are negligible. Also the overall 
structure is slender. So the stiffness of the overall tower 
both for space truss and frame model remain close. Hence 
as a first step, the tower is modelled as a space truss for 
the eigenvalue analysis. 

The tower studied has 6 p^els in the body and the 
base width of the tower is taken as 4,33m, shown in Pig. 2,1, 
Lumped mass formulation is used. Pig, 2,2. Pirst 6 natural 
frequencies (square- root of eigenvalues) of the tower are 





Fig.2-2 Tower idealization 
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extracted using subspace iteration method. The results 
are presented in Table 2,1, The first eigenmode is the 
bending mode in the transverse direction, second is the 
bending mode in the longitudinal direction of the line and 
the third is the torsional mode. The fourth, fifth and 
sixth natural frequencies are very well separated from the 
first three. This shows that the major portion of response 
mainly depends on the first three modes , 

T abl e 2 , 1 ITatu ral gr® quen cj.es^ pf ^Tr^jmis^sJLro Towe r with _ 
15m Body Height 


Mode Ho, 1 2 3 456 


Structure 

frequency 

(rad/sec) 26.3 26.6 35.1 77.7 84.0 89.6 


Coming to the individual members, the long member 
is the diagonal in the bottom most panel . The length of this 
member is about 5.25m. Considering this member as simply 
supported, the lowest natural frequency is given by 



P A 


( 2 . 1 ) 


where m = 


g 
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Substituting the values, P= 0,00783 kgf/cm^, 

2 2 

E = 2,04-7 X 10 kgf/cm , g = 981 cm/sec and the maximum 
value^^ of (1/r) = 200, we get 

*^.1 = 4-7.75 rad/sec 

Since this diagonal is subdivided by the other diagonal 
in the same panel (also by secondary braces), on the 
face of the tower, the lowest natural frequency become 
w-j = 4 X 47.75 = 191 rad/sec. But in the perpendicular 
plane to the face of the tower, however, the member 
frequency is only 47.75 rad/sec. Hence the lowest natural 
frequency of the long member is well above the first three 
natural frequencies of the tower, Therefore, the response 
of the overall structure shall govern the design rather 
than the response of individual members. 

It is concluded from the above discussion that the 
transmission tower need be modelled as a space truss, GHie 
space truss model result in the number of degrees of freedom 
same as that in static analysis . 

2*5 Oonfiguration 

2,3,1 Basic Bracing Configuration : 

There are three types of bracing systems which have 
been mainly used on tower bodies (Fig, 2', 3) , 
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( a) Tension system, 

(b) X- braced tension-compression system, and 

(c) K- braced tension-compression system. 

In the tension bracing, all the diagonal members 
have slenderness ratio high enough to act only in tension. 
Such bracing results in larfee deflection under heavy loading, 
because the tension members are smaller in cross section 
than v^at the coimpression members would be for the same 
loading. Overloading of such towers would result into 
the failure of the inactive tension member in compression, 
and the active tension members would carry the load and 
prevent the collapse of the tower. X-type of bracing 
arrangement is very efficient for relatively small lateral 
dimension of the tower. Such arrangement gives better 
distribution of load to the tower footing, K-type is 
generally applicable to large structures and has not been 
found competitive with either the tension system or the 
X-bracing for normal suspension towers. 

It is possible, in general, to combin© two or more 
of the above basic types to obtain a combined bracing 
configuration for a particular tower, which may prove to 
be economical'. However, for the purpose of this study, it 
has been assumed that all the panels have X-type of bracing 
which has been found to be economical for the body heists 
upto 25 m^^ . 
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2 ,3 .2 G-e neration of Oonfi ff Tatlon 

The confi.'mration. of the tower is generated In two 
parts. The first part generates the fixed data such as 
connectivity of the memhers and t 5 rpe of members, which are 
computed only once during the optimization cycle. The 
second part generates the coordinates of the joints and 
lengths of secondary members which are generated after 
any change in the design variables. 

2, 3. 2.1 Generation of Fixed Data: Subroutine ’DATAHT’ 

The program has been written in general so that 
data for any latticed tower structure can be automatically 
generated with the minimum possible information as input. 

The tower structure has joints interconnected by members, 
located on horizontal planes at various levels . At such 
a plane lies four or eight joints as shown in Fig. 2,4 
depending upon the type of bracing adopted, i.e., X-type 
or tension bracing result in 4 joints in a horizontal plane 
while K- braced panels result into 8 joints. These joints 
are located equidistant from the vertical axis of the tower 
and hence if the coordinates of one joint (henceforth called 
as reference joint) is Tmown, then coordinates of the other 3, 
in the case of X-bracing, can be generated. The joints are 
nximbered sequentially in the anticlockwise direct ion, starting 
from left joint on the front face of the tower, in each level. 
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I'or exp’nple, the coordinates of (1) in Fig. 2,4 will only 
he fed and the coordinates of (2), (3) and (4) will he 
generated dimply by changing the sign of the coordinates 
of (1) * If a panel has K- bracing, then the coordinates 
of the middle Jaint on the front face of the tower is 
also required to he fed, 

Ihe other information required for the tower 
configuration is the connectivity table of the members. 

Far this purpose the members of the tower are grouped as 
leg members, diagonal members, horizontals etc., in each 
panel, lo generate the connectivity table, the connectivity 
of only one member from each group is read and for all the 
other members, the connectivity is generated based on 
symmetry of the tower. This is pos-sible by suitably 
numbering the joints as shown in Fig, 2,4, For example, 
connectivity of leg member 1 will be fed as data while the 
same for other members 2, 3 and 4 is generated by suitably 
incrementing the joint numbers of member 1 . Similarly for 
diagonal bracings, connectivity for member 5 is fed while 
the data for other 7 members in this group are generated. 

In Fig, 2,4, the members are numbered starting from 1 for 
the sake of illustration. The computer program written 
is, however, general and all what is required is that the 
numbering of joints in the same sequencial order. 



34 


Por K- braced panelB, the member data is generated 
in the similar way. The only difference is the increase 
in the nmber of groups in a panel. It is 4 for K- braced 
instead of 2 for braced panels. The typical input data 
for the panel shown in Pig. 2,4 is given in Table 2.2, 

Table 2 ,2 Typical Data for Q-enerating the Configuration 


X-bracing 

K«- bracing 

! 

J 

member 

number 

IT* 

IP(J) 

IQ(J) 

A(J) 

j 

IT 

IP(J) 

IQ(J) 

A(J) 

1 

1 

1 

5 

# • « 

1 

1 

1 

9 

• • 

5 

3 

5 

2 


5 

2 

9 

5 

• # • 





i 

13 

5 

5 

2 

• « • 
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4 

5 

6 

• « • 


* IT = 1 for leg members 

- 2 for K- braces 

= 3 for X- braces 

= 4 for horizontals (4 members) 

in 

= 5 for horizontals is K-braced panel (8 members) 
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The coordinates of the joints in the hashet portion 
of the tower is fixed and hence generated only once. The 
connectivity table of the members of the tower remain the 
same, except that the areas of the member goes on changing. 
When the area of a member is changed, the latest value will 
be stored in place of the old value. 

Other data read in this subroutine are the width of 
basket, height of tower body, and height of panels and base 
width of starting design point, 

2. 3.2 ,2 G-eneration of Variable Data: Subroutine 'XIZCOR* 

In the present work, the configuration of the 
tower body changes during the process of optimization since 
the panel heights and base width are the design variables. 
Hence the coordinates of the reference joints changes from 
design to design. To start with, the coordinates of the 
reference joints are known which depend upon the starting 
point of the optimization subroutine. In otherwords 
the coordinates of the reference joints can be worked out 
once a set of values for panel heights and base width is 
known, Knowing the width of basket, width of base and the 
height of tower body, the slope of leg members in the body 
can be worked out. With this slope and height of panels, 
the coordinates of the reference joint can be computed. 

Once the coordinates of the reference joint is known, the 
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Coordinates of the other joints in that panel can he generated 
as discussed in the preceding section. 

Ifter generating the coordinates of all the joints 
in the body, the program computes the length of secondary 
members in each panel. If the length of secondary brace is 
less than 2,17m, thinnest possible prescribed section is 
adopted. The area of secondary braces, whose length exceeds 
2,17m, is piched up from the list of equal angles so that 
slenderness ratio of the section taken is less than 250, 

The weight of all these secondary braces is summed up and 
is added to the weight of primary members to get the total 
weight of tower at any design point, 

2,4 Static loads Acting on the Tower 

2,4,1 . S .Codes f IS: 802 and IS; 875 

The transmission line towers are situated mostly 
in open country. Hence one of the main loads on the 
transmission line towers is due to natural wind. The natural 
wind force depends on many factors which include local 
topography, climate variation, height above ground etc , and 
as such the wind load should be treated under probabilistic 
framework. The importance of this necessity can be understood 

I 

from the recent works on wind loading of structures reported 
in literature^®“^^ . 
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To treat the wind load under prohabilistic framework, 

it is necessary to have statistical data of wind pressures 

observed around the site under consideration for a considerable 

period of time. In India, such data is not available, but 

66 

instead Indian Standard, IS: 875-1964, has recommended for 
general structural design wind pressure values for different 
locations of the country* India is divided into three 
zones for the purposes of presenting the wind map of India. 

The wind pressure values, including winds of short duration 
as in squalls, are presented in the form of table for 
increasing hei^t. The basic wind pressure indicated is 
the ma:simum ever likely to occur in the respective areas, 
under fully exposed conditions . 

The Structural liigineering Sectional Committee, 

constituted in 1967 to draft the code on »Use of structural 

steel in overhead transmission line towers: Part I- Loads 
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and permissible stresses' has recommended wind pressure 
values for the design of transmission line towers -based on 
IS: 875-1964, The wind pressure map and the table giving 
the pressure values, presented in IS: 802-1967, were taken 
from IS: 875* Guided by this code on transmission line 
towers, the towers in India were designed till mid- seventies , 
This code has prescribed different factors of safety for 
the design of structural members under normal condition in 
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which all wires are intact and broken wire condition during 
which one of the conductors or ground wire is assumed to be 
broken* The details are presented in the next section. 

later, the sectional committee felt the need to 
introduce the dyiiamic effect since wind induce considerable 
vibration in the tower due to the variation of wind pressure 
with time. This points to the fact that the wind pressure 
values to be taken are static equivalent of dynamic wind 
rather than maximum static pressure; mflaen the code IS: 802 
(part I) was revised in 1973 the committee has recommended 
an increase of about 30 percent in wind pressure values 
that are presented in IS: 802 (part I)'-1967, irrespective 
of height of the tower. This does not seem to be logical. 
The stiffness of the tower varies with height and plays 
an important role to decide the static equivalent of 
dynamic wind pressure. In other words the increase in 
static maximum wind pressure, to account for the dynamic 
effect, has to be arrived at for different heights of 
tower and cannot be a constant. Hence the design of 
transmission line tower taking wind pressure values from 
IS: 802-1973 would result into very conservative design 
especially for smaller towers. This makes it necessary 
to analyse the transmission line tower under dynamic wind 
pressure. Only the dynamic analysis which takes the effect 
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of stiffness and mass of the tower into account, and wind 
as a function of time would produce rational design. 

Based on the above discussion, loads for the static 
and dynamic analysis of the tower is derived from 13:802-1967. 
Comparison of results of Chapters 4- and 5 gives the effect 
of d 3 niamic wind over the static loading. 

The latest code, i.e., IS: 802-1973, refer the 
intensity of wind pressure in the three zones as light, 
medium and heavy. The factors of safety recommended remain 
the same as in the former code, hut there are changes in the 
allowable stresses and allowable maximum slenderness ratios 
for the members of the tower* Calculation of tower loadings 
for a typical 132 K7 double .circuit line tower is also 
presented in the latest code, 

2 ,4 ,2 Computation of Loads IS : 802 

Transmission line tower is subjected tb three types 
of loads^®^^^ viz,, (1) vertical loads, (2) transverse loads, 
and (3) longitudinal loads. The vertical loads are computed 
on the basis of weight span, shown in Fig. 2.5, which is 
the distance between the lowest points in the catenary in 
the spans adjacent to the tower. Transverse loads are due to 
wind on bare or ice coated wire and also due to change in 
line deviation. The computation of wind load is based on 
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The loading, on the 220 KY double circuit tower (Pig. 1.1), 
for different cases of normal and broken wire conditions 
presented in Pig. 2.6 which are taken for the present work. 
These loadings are obtained from Cromptons India Ltd. The 
wind direction is reversible and hence the total number of 
independent loading conditions turn out to be 8 as follows: 

1. Normal condition, i.e., all wires intact and wind 
blowing. 

2. Broken wire condition of ground -wire and wind blowing. 

3. Broken wire condition of top conductor and wind 
blowing, 

4 . Broken wire condition of middle conductor and 
wind blowing, 

5 * ; Borkenjwire condition of bottom conductor and c;! — 

wind blowing. 

6. Same as case 5 but wind in opposite direction. 

7 . Same as case 4 but wind in opposite direction, 

8. Same as case 5 but wind in opposite direction. 

However, a simple study reveal that the number of critical 
loading conditions are much less. This reduces the time 
required for analysis and hence is a substantial gain for 
carrying out the optimum design, 

2,4,3 Critical Lo a d Conditions 

Por the study of critical loading conditions, a 
tower with 15m body height having 63 joints and 208 members, 
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(Broken wire loads are given in parentheses) 
( Loads in kg & dimensions in m ) 


Fig.2'6 Loads on the tower basket 
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shorn in Eig. 2.7, is analysed subject to all the eight 
loading conditions as discussed in the preceding section 
and the self weight as the ninth case, through this study 
the following observations are made, 

( i) Self weight has considerable effect on lower leg 
members. This increases the compressive force in the 
lower leg members on the leeward side and hence has to be 
considered for analysis to be embedded for the automated 
optimum design. This is due to the fact that automated 
optimum design seeking a minimum weight tower endeavours 
to pick up member areas which will lead the constraints 
to be bounded. 

(ii) Normal condition need not be considered, as this 
does not produce worst condition of loading for members 
except the horizontal members at the junction of the 

cross arms . Ttirthermore , maximm stress in these horizontal 
members are at a level for which minimum section prescribed 
for tower members are safe, 

( iii) The broken wire condition of bottom conductor is 

not critical for any member except those two group of members 
meeting to foim this cross arm. It is, therefore, not 
necessary to consider this case, for the purposes of opbimum 
design of tower, since the increase in weight due to broken 
wire condition of this conductor shall be a constant , 
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( iv) Some members (bencefortb termed a group in this 
section) are critically stressed in compression when the 
wind is blowing from left to right as shown in Pig.2.6 
and one of the conductors on the left face snaps. It is 
observed that when the direction of wind is reversed, then 
the members in the group experience maximum tensile force 
whose absolute value is less than the maximum compressive 
force experienced in the earlier case. Moreover the 
maximum compressive force in the group in the latter case 
is less than the maximum compressive force experienced 
in the former case. Therefore, it is good enough to 
consider broken wire condition with only wind blowing 
frcan it towards the vertical centre line of the tower. 
Member forces produced due to the eight different loading 
conditions listed earlier are tabulated for typical members 
in Table 2.3 in which the values that govern the design 
are underlined. The member forces listed in Table 2.3 
are obtained by multiplying the internal force computed , 
through analysis by the corresponding factor of safety 
for the loading condition considered. Por a few diagonal 
members the critical compressive force is produced due to 
loading number 7 or 8. In most of the cases of diagonal 
braces, this does not govern the design. It is the l/r 
ratio restrictions, where 1 is the effective length of the 
member and r is the least radius of gyration, which govern 
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in the same order as listed in Section 2,4.2 
compression , + tension 
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the' memher size . 

Based T>n the above observations it is concluded 
to analyse the tower only for the three cases of loading 
listed below, 

(1) Ground wire broken condition and wind blowing. 

(2) Top conductor broken and wind blowing from it 
towards the vertical centre line of the tower. 

(3) Middle conductor broken and wind blowing from it 
towards the vertical centre line of the tower. 

The bajsket portion of the tower with the above mentioned 
critical loading conditions, in addition to the loads 
corres^nding to the normal condition, is shown in i'ig,2.6. 

2 .5 ly ynamic loads Acting on the Tower 

\loads on the structure are, in general, dynamic, 
with an exception of dead loads. Hence the idealization 
of structural loading, as a function of time is an important 
step in the present work. 

The loads on the transmission line tower can be 
classified into three groups, viz.,(1) dead loads, (2) broken 
wire loads, and (3) wind loads. The dead loads, of course, 
are static loads'. The conductor or ground wire breaks 
during a turbulant wind, say, during heavy cyclonic storm. 
Hence these broken wire loads are dynamic loads,. 'What 
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precisely is the load on the tower 'in this situation is 
rather complicated^"^. However, IS: 802-1973 has prescribed 
certain guidelines to treat these broken wire loads as 
static loads. Hence the broken wire loads are treated 
like static loads in this work as described in 
Section 2.4,2, The wind load acting on the tower is a 
time dependent forcing function and its idealization is 
the subject of discussion in the following section. 

2 ,5 ,1 ' Wind Loads 

The wind that passes over the transmission tower 
is a lange scale movement of air, due to thermal currents 
and the rotation of the earth. The wind velocity which 
has temporal as well as spatical randomness, depends upon 
many factors like pressure, temperature, the particles 
carried by wind etc. In order to design the structure 
against wind loads a continuous record of temporal and 
spatial variation of wind, at the site of construction, is 
essential. But more often than not, this is not possible 
because of a multitude of reasons. In such a situation, 
some alternative approach has to be adopted so as to avoid 
the need of having a continuous wind records. 

68 

Maharshi has adopted ‘Digital- simulation* technique, 

in the absence of wind records, to analyse a 175m high T.V. 

6Q 

tower. The simulation has 'been done with Chiu’s ^ model as 
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- 70 

well as with Shinozouka” s model. His work was mainly 
to study the effect of parameters with different mean wind 
velocities and compare with that of static wind loading 
taken at mean value. 

To carry out the dynamic analysis with simulated 

(pseudo-wind) records, the required informations are, 

(1) mean wind velocity, (2) the standard deviation of 

wind velocity, and (3) time interval during which the 

wind is assimed to vary linearly. Even these are not 

71 

available for most of the locations where failure due 
to cyclones are estimated heavy. Originally the author 
was interested to study the response of transmission tower 
to that cyclonic storm which hit the Indian east coast 
during November, 1977, resulting in huge loss of life 
and materials. Quite a number of transmission line towers 
collapsed in that cyclone. The Indian Meteorological 
Department, when requested, could provide only the peak 
values of wind speed observed which are given below. The 
maximum wind speeds, given below, correspond to the Bay of 
Bengal storm of 14.11 .1977 to 19.11.1977. 

(1) Maximum wind speed reported by ship on 17.11.1977 is 
about 90 knots, 

(2) Maximum wind speed estimated from satellite pictures 
is about 120 knots on 18th and 19th of November, 1977. 
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(3) Maximum -wind speed recorded at Gannavaram (Andra 

Pradesh) when the storm moved close to the station 
on 19 . 11.1977 is 75 knots. 

They expressed their inability to provide any wind record 
or the statistics of wind speed variation. The above 
informations are as good as that available in IS: 802-1967, 
since no information regarding the variation of wind speed 
is available. 


With the limitations discussed above, it is decided 
to Carryout the deterministic analysis, in which the wind 
load is approximated to a known function. 

72 

Yeh and Shieh have represented the wind load, 

for the analysis of hyberbolic cooling towers, by a 

deterministic load function. The load function assumed 

has the peak values, at certain time intervals, representing 

the gust. The load-time history used in their work is shown 

in Pig. 2,8. The numerical value of peak wind pressure, 

taken by them, is of no importance here and hence the peak 

73 

value is represented by I'^iter Steinmetz et al have 

studied the behaviour of hyberbolic cooling tower for which 
wind loading was based on full scale measurements by pressure- 
difference transducers. The pressure readings were taken 


at the throat of an existing tower located in Martinis creek, 


Pennsylvania. The d 3 niamic pressure has 


I*' 4- '.r 
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Courier series. They found, on ooi^arison, the basic 
response patterns were similar to Yeh’s tower. This shows 
that loadv-time history shown in Fig, 2,8 can be used to 
represent wind without much of error. Hence the deterministic 
loa<d-time history shown in Fig, 2,8 is adopted for the 
present study. 

Wind pressure varies with height and hence the 
maximum wind pressure value P^goj; in. Fig, 2,8 depends on the 
height above ground. Based on the wind pressure values in 
IS: 802-1967, the maximum ■ dynamic wind load, for the . 

t 

! 

transmission tower analysis, has been computed. However, I 

;6kot«n I 

the time history of wind, s te am in Fig, 2,8, is taken as | 

same for all heights above ground, i,e., time interval of > 

I 

peaks and the level of mean value of wind load remain the | 

same for all dynamic wind loads that are assumed to act at j 

i 

joints of the tower. The wind loads, computed at maximum I 

wind pressure from IS: 802-1967 as well as the other | 

static loads are shown in Fig, 2,6 , In this figure all I 

i 

the loads in the X-direction (wind direction) are dynamic 
loads for which the time history is shown in Fig. 2.8. 

i 

2 ,5 ,2 Esirthquake loads 

i 

IS; 802 has not mentioned anything about the j 

consideration of earthquake in the design of transmission 
tower. However IS: 1893-1975'^^ has given certain guidelines f 
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for the design of structures in general . The average 
acceleration spectra given in the above code, is taken for 
the present work of transmission line tower design. The 
code has stated the following assumptions applicable for 
general structural design. 

(1) Earthquake causes impulsive ground motion which is 
complex and irregular in character, changing in 
period and amplitude and lasting for small duration. 
Therefore, resonance of the type as visualized under 
steady state sinusoidal excitations will not occur 
as it would need time to build up such ampli Ludes , 

(2) Earthquake is not likely to occur simultaneously 
with wind or maximum flood or maximum sea waves , 

The country is classified into five zones for the 
purpose of determining the seismic forces; the details are 
available in IS: 1893-1975. for the present study the tower 
is assumed situated in the most critical zone, i.e., Himalayan 
region, zone 7, 

Depending on the structure, the code has recommended 
two methods for computing the seismic force: 

(a) Seismic coefficient method, and 

(b) Response spectrum method. 

In the present work, the response spectrum method is used, 
which takes into account the system properties explicitly. 
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where N is the mass and g is acceleration due to gravity, 
This horizontal force is assumed to act at the corresponding 
location of the mass. For the present study the direction 
of seismic acceleration is assumed in the transverse 
direction of line. Further details of modal analysis and 
dynamic stresses in the members of the tower due to 
earthquake are presented in Chapter 5 . 





OHiPTER 3 


-FOmULAUllOE OP OPTIimi DESIGN PROBLIM AND METHOD OP OPTIMIZATION 

5 • Intro duc t i on 

’ii/hen a means for predicting the behaviour of any 
structure is available, and design philosophy is decided, 
it is possible to cast the design modification problem 
in the form of a mathematical programming problem, 

A mathematical programming problem is one in which 

a multivariable function f(d) , where d is a n ^dimensional 

vector consisting of d., 3 = 1 , 2 ,.., n, is minimized (or 

J 

maximized) subject to given constraints gj_(d) 

i=1,2,,,,,m (m _< n) , where b^ eire known constants. There 

are usually inequality constraints on some or all of d.. 

J 

The function f(^) is called the objective function and its 
choice is governed by the nature of the problem. Por a 
configuration optimization problem, cross sectional variables 
of member as well as the variables that define the geometry 
of the structure are included in vector d. 

The objective of stiuctural design optimization is 
frequently taken to be the weight minimization. This is 
justified for transmission line towers which consist of 
sections made of only one material (generally steel) as well 
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as the fact that fabrication and erection cost is on 
tonnage basis. The minimum weight design is, therefore, 
considered as objective function in the present work, 

5*2 Formulation of the Optimization Problem 

The general structural optimization problem can be 


stated as follows: 



Minimize f(d) , d^ 

= ^2 V 

(3.1) 

Subject to 



iY(d,x,t)| < 

y(u) 

^ t > 0 

(3.2) 

la(d,x,t)| < 



< V < 

i ^ 

u) ^ _ -j ^ 2 , , . , , p 

(3.3) 


k = 1,2, q 

(3.4) 

and < '1-i < 

d^^^ j = 1,2,.,,,,m 

J 

(3.5) 


where f represents the weight of the structure 

d vector of design variables, viz,, geometrical 
dimensions and cross sectional variables 
I Y(d,x,t) I represents the displacement at any point on the 
structure as a function of space variables x and 
time t 

y^^^is the upper limit on the displacement 
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la(d,x,t)l 


0 


(u) 


xl'^) 

1 

and 


represents the stress at any point on the 
structure as a function of space variable x 
and time t 

is the upper limit on the stress, taken to he 
the maximum permissible stress as per design 
standards 

represents the i-th eigenvalue of the structure 
is the upper limit on the i-th eigenvalue 
is the lower limit on the i-th eigenvalue 

is a function of configuration variables only 

and hence d is a subset of vector d 

constant, known from the compatibility requirement. 


The behaviour quantities Y, a , and are dependent 
on the design variables and hence the constraint Eqs.(3.2) 
and (5.3) are called behaviour constraints. The compatibility 
of structural configuration is imposed by Eq. (3.4). The 
side constraints, Eq. (3.5), impose the limits on the size 
of the design variables based on minimum size of sections 
and any restriction on the geometry of the structure. 


For the static loads, the time parameter is set 
to zero and hence the problem gets simplified as follows: 


Minimize f(d) , ^ = Id^ d 2 



• • • « 
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Subject to 

Y(d,x) < 

o(d,m) i 

(3.6) 


...,q 

(3.7) 

and 

i 3 = 1 , 2, ....,m 

(3.8) 


where vaxious quaiitities have same physical meaning as 
described earlier. 


5.2.1 Design Variables 

As stated in Chapter 2, the basket portion of the 
transmission tower geometry is defined from electrical 
considerations. Hence this portion of the tower geometry 
has little scope for variation with the aim of reduction 
in weight, 

The geometry below the bottom cross arm, ie., 
body of the tower, is free from any restriction and hence 
it has ample scope for modification. All that is required 
is the total height of body should be kept constant through 
out the process of optimization. The configuration of the 
tower body is defined from the width of base, number of 
panels and height of each panel. Hence the design variables 
chosen in the present work are the panel heights of the 
tower body and the base width of the tower i J'ig. 3.1. 
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Available angle sections are used for the members 
of the tower. Hence, for the design of members the req.uired 
section is picked up from the list of available sections 
and, therefore, the cross sectional variables do not appear 
explicitly in the optimization process. The details are 
discussed in Chapter 4. 

5.2.2 Objective Hunotion 

The goal of the present work is to study the effect 
of panel heights and the base width on the minimum weight 
of the tower. Hence the objective function, f, representing 
the total weight of tower, which is to be minimized, may be 
stated as follows: 

V 

f =P 2 A. 1. n. , i = 1 to n (3.9) 

X i X g 

where P = density of the material (steel) 

1^ = length of members in group i (group defined 
earlier in page 33 ) 

= area of members in group i 
= number of members in group i 
and n~ = number of groups which is the sum of groups 
of primary as well as secondary members, 
l^s are length of various members and are implicit function 
of design variables which are: 
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(1) tlie panel height h^, i=1 to (KP-1), NP being the 
number of panels in the tower body as shown in 
Pig. 3.1. Only (ITP-I) out of NP panel heights are 
treated as explicit design variables to make the 
formulation unconstrained. The details are discussed 
in the next section, 

(2) the base width of the tower, b. 

A^s are also implicit functions of the design variables, 
since they also depend upon the configuration of thd tower. 

3.2,3 Constraints 

The constraints that are to be satisfied, in a 
general optimization problem, have been listed earlier. 
However, the following assumptions in the present formulation, 
under static and dynamic loading conditions, make the 
optimum design problem unconstrained, 

1, The displacement and stresses in the structure are 
very well correlated, i,e., any constraint on stresses 
lead to imposition of constraint on displacement implicitly. 
Further, constraint on displacement is not an important 
factor in transmission line towers so long as it does not 
undergo excessive deflection, due to instability of the 
tower. Therefore, restriction on deflection is not 
necessary. Hence, only constraint on stresses have been 
considered in the present work. The manner of handling 
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these constraints, implicitly, is discussed in Section 4.3. 

2, In general, the constraint on natural frequency 

is imposed so that the optimum design obtained has natural 

75 

frequency away from resonant natural frequencies . However, 
any design that falls in resonant region will require more 
weight. Even if the formulation does not mpose any 
restriction on natural frequency, the program will move 
away from the resonant region, since it is searching for 
a minimum wei^t. The results in Chapter 5 brings out 
this aspect very clearly. Hence in this formulation no 
restriction on natural frequency is imposed, 

3, The body height of the tower is always constant, 

Bq. (3.4) and Eq. (3.7); i.e., we are attempting the 
optimum configuration for a given body height of tower. 

This equality constraint can be stated as, 

h^ + h 2 + ~ constant (3.10) 

where h^^ is the height of i-th panel from top in body. 

To handle this constraint implicitly, only the first (HP-I) 
panel heights are taken as design variables in addition 
to base width, making the bottom most panel dependent on 
the value of other panels, 

4, Since the available angle sections are used, the 
section having minimum as well as maximum area stand for 
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the hounds imposed hy Bq, (3.5)* The details are discussed 
in Section 4.3. Hence all the restrictions on the member 
selection, stipulated by IS: 802>-1973, are handled implicitly, 

3.2,4 Mathematical Hature of Bormulated Proble m 

The optimization problem formulated in the previous 
section can be stated as follows: 

Minimize f(d) , where d = [h^ hg .... "b ] » 

since all constraints imposed are being handled implicitly. 

The unconstrained function, f, which is the weight of the 
tower to be minimized is not an explicit function of the 
design vector, d . However, for a knovm. design vector 
function value f can be calculated by carrying out the 
analysis and the design of the tower. Thus the formulation 
of the optimum configuration design of transmission line 
tower turns out to be an unconstrained minimization problem. 

5.3 Method of Op timization 

Several methods are available for solving an 

unconstrained minimization- problem . In the case of simple 

functions one can discuss the efficiency of a particular 

method over the other. But for the engineering applications 

the nature of problem makes the choice limited. The paper 

77 

by Carpenter and Bmith gives a comparative study of 
unconstrained minimization algorithms, using in the Sequencial 
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7R 

Unconstrvained Minimization Technique (SUI4T) , to optimize 
a variety of structural engineering prohlems. While listing 
their conclusions they have rightly pointed out that 
conclusions concerning computational efficiency of optimization 
algorithms can never he stated categorically because an 
exceptional problem can invalidate conclusions drawn from 
other investigations. 

The unconstrained minimization methods can be 
classified into two broad categories as non- gradient methods 
and gradient methods. The non- gradient methods require 
only the evaluation function values. The gradient 
methods require, in addition to function values, the 
evaluation of first and possibly higher order derivatives 
of the objective function. 

The evaluation of gradient could be done using 
finite difference technique, since the function to be 
minimized is not an explicit function of design variables. 
However, this entails additional function evaluations. 

J^art from the time consumed for such function evaluations, 
since the area of angle sections that are used for the 
members of the tower is discrete in nature, the gradient 
approximation may not be accurate. Hence non-gradient method 
is the best choice for the problem foimulated in the present 
work. 
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Among the non- grad lent methods, Powell’s method ^ is 
the most widely used search method since it is quadrat ically 
convergent, Therefore Powell's method is chosen in the 
present work for the unconstrained minimization, 

3.3.1 P owell's Method 

The sequential steps involved in the Powell's method 
is presented in the form of a flow chart in Pig. 3.2 and 
discussed briefly below. 

(1) Take a starting point and linearly independent 
search directions. To start with, the search directions are 
taken to be coordinate directions and hence are linearly 
independent. 

(2) Find the step length a* such that the objective 

function has minimtim value, i.e., ®n’ 

searching along n-th coordinate direction. This design 
point x-j is stored for later use. 

(3) For i = 1,2, ..., iterate the following steps(4)-(7) , 
till convergence is achieved. 

(4) For r= 1,2, n, find, along each direction, 

such that f(^ + oc* s^ is minimum and update the design 
^ 

vector as x^^ ~ ^ “r 

(5) Generate the pattern direction given by s^^^ = . 

The design point x^^ is the latest point in step (4). 
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(6) Pind oc^ such that Sp ) is minimum, 

store ^ + “p Sp as x^. Since is a iDetter direction 

than coordinate directions discard one of the coc^rdinate 
directions, For example in the first cycle would he 
discjurded. Renumber the direction vectors, as the first 
direction is discarded, as follows. 


s 


i 


’i+1 


and Sj^ 



(7) Check for convergence. If optimum is reached, come 
out of the cycle. 


In the above procedure, the routine moves through 

s^ direction twice in the first cycle, steps(2) and(4) , 

and obtains design points x^ ^d Xj^ which are minima along 

s^ and hence the direction s^**^ and similarly subsequent 

pattern directions have special properties that are the 

source of rapid convergence of the method. The details 

79 

are given in the Powell’s paper as well as in books on 

... .. 76,80 

optimization * • 

In step(7)of the procedure, the routine checks 
whether optimum is reached, Powell recommends a termination 
criterion for ordinary use such that when a cycle produces 
a change in all variables of less than one-tenth of the 
required accuracy, the process is stopped, A. safer (i.e., 
less likely to stop prematurely) criterion is also given by 
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79 

Powell which is much more time consuming. 

In the Powell's method, it is necessary to find a 

optimum step size a in each search direction (step 4)* 

For this linear minimization, quadratic inteipolation 
80 

technique is used in the present work as this method 
does not require any information about the derivative of 
the function. 

In the present work, to find a useable search 
direction, a step size of S = 0.01 has been used, to 
check whether the function decreases along this dicection, 
i.e., f(x + 3s) < f(x). If it is not a useable direction 
with positive value of 6, 3 is taken as -0.01 and the process 
is repeated. It has been observed, in certain situations, 
that the function values at x + 6s as well as at x - 3s 
are found to be greater than f(x). However there is 
difference among the function values evaluated on both 
sides of the present design point x , This information is 
used to establish the sea^^^ch direction, even though it is 
a violation of basic assumption that the objective function 
is unimodal. Since the available areas of the angles are 
discrete in nature, this type of search has been attempted 
and found useful in few directions. 

In the search along a useable direction, objective 
function values are evaluated at a =0,5, 0.6, 0.9t to 
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bracket tbe minimum ; i.e., if tbe direction of search is 
the coordinate direction, then changes in design variables 
are 0.3 m , 0.6m,.., etc. Of course, if f(x^ + 0,3 >f(^) 

then a = ^ etc,, is attempted, ^len 

a becomes less than 0*02, in the above process, then 
further search in this direction is terminated and the 
last value taken as the minimum along this direction, 

The optimization process is terminated when the 
change in all design variables is less than 0,001 in the 
search along the pattern direction. 

The first (n-1) variables of the design vector 
correspond to panel heights starting from top and the n-th 
design variable correspond to the base width (Pig. 3.1). 

Hence to start with , in the first cycle of Powell’s method, 
the routine finds a better base width locally, goes on to 
panel heights and finally moves along the pattern direction. 



CHAPTER 4 


0Ptl!([U5i[ DBS ten POR STATIC LOADS 

^ jlntj c o dgctiQn 

Static loads acting on the tower structure have 
been discussed in Chapter 2, In this chapter, the analysis 
of tower as an elastic space truss subjected to those 
static loads is described first. The discrete element 
method of structural analysis is used in the present work. 
The discrete elements refer to the primary members of the 
tower only. The secondary bracing members of the tower 
are not considered for analysis. However, the effect of 
secondary braces which subdivide the length of primary 
members are taken into account for the desigii of primary 
members. 

The optimum configuration design under static 
loads is carried out using Powell’s method described in 
the previous chapter. The resulting optimum dimensions 
for different heights of tower body and the discussion 
on the results are presented in this chapter. 

4- • 2 Methods of Static Analysis 

The analysis of the tower subjected to multiple 
loading conditions .form a significant part of the total 
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effort required to design the tower. Optimum design process 
being a recursive analysis- design procedure, one has to 
perform a large number of analysis of the tower before 
arriving at the optimum design. Therefore, the choice 
of an efficient method of analysis is warranted. Hence 
a study has been made in the present work to choose from 
among the flexibility and the stiffness methods of structural 
analysis. 

^ flexibility Method 

A method of analysis useable in design office is 
given in Ref, 81 which presents a systematic procedure 
for establishing accurate distribution of shear forces, 
arising from torsional moments, on to the faces of the 
space tower. However, this procedure, which falls in 
the broad category of flexibility method of analysis, is 
considered not suitable for computerization in the present 
context. 0?he flexibility method that is programmed for 
the present purpose of automated analysis is presented 
below. 

The basic procedure for analysing a structural 
system by the flexibility method is as follows; 

i) Release the original structure so that it is a 

stable determinate structure; the released actions 
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(reactions and/or internal member actions) are 

called redundants, 3 “ 1 , 2 ,.... n , where n 

3 r 1 

is number of redundants. 

ii) Calculate the displacements, at the locations of 
released actions, due to the original loads in 

the direction of the released actions (redundants). 
Let these displacements be put in vector form d^^ 
where subscripts q and 1 stand for redundants and 
applied loads respectively. 

iii) Jipply unit values of the redundants, one at a 
time, to the released structure. 

iv) Calculate the deformations at each release in the 
direction of the release due to the unit value of 
each redundant acting separately. 

v) ITite a compatibility equation for each released 
action by describing the displacement of the point 
of application of the redundants as a linear 
combination of the deformations calculated in 
steps (ii) and (iv), 

vi) Solve the set of linear algebraic equations for 
the unknown redundants. 

vii) Place the known values of the redundants back on 
the stiructure and complete the solution by using 
the equations of static equilibrium. 
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The equation of compatibility, in matrix form, 
along the line of redundants can be written as 

^ 1 + ^ql = 0 (4*1) 

where P is the flexibility matrix and q is the vector 

of redundant forces. Eq, (4.1) has S3miultaneous 

equations corresponding "to redundants, Eor multiple 

loading conditions, the vectors q and d^^ become matrices 

of size n^ x n^, where n^ is the number of loading conditions. 

Details of this method can be obtained from standard boohs 

82 

on structural analysis . 

4*2.2 S tiffness Meth od 

The stiffness method is being used extensively to 
solve structural analysis problems due to the elegant nature 
of the method from programming point of view. Basically 
this method consists of writing the equilibrium equations^^ 
at all the joints in terms of the displacements as 

K u = p (4.2) 

in which K is the assembled stiffness matrix of the structure, 
u is the vector of unhnown displacements at the joints and 
p is the knomi vector of external forces acting on the 
structure. The displacements can be obtained by inverting 
the Z matrix and post multiplying it by p. Moreoften, 

Eq, (4.2) is solved as a set of simultaneous equations®'^ 
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for fhe unknoms, j = 1,2, where is the 

number of degrees of freedom. Solution as simultaneous 
eq_uations is less time consuming due to banded- symmetric 
nature of K matrix, Knowing the vector u , the internal 
forces Can be obtained from element stiffness matrix®^ and 
from the components of u, For multiple loading conditions, 
vectors u and p become matrices of size n^ x n^* 

4 • 2 , 5 Results leading to Choice of Method of Static Analysi s 

The tower shown in Fig, 2,7 has been considered 
for the study of suitability of flexibility or Stiffness 
method viz- ar- viz computer storage and time requirements. 

It has 31 degrees of indeterminacy. Incidentally all these 
redundants fall in the basket portion of the tower. Fig. 4,1. 
The body of the tower has no redundant member when Z-type 
of bracing is adopted. 

The program starts from ^oint 1 and goes on solving 
the static equilibrium equations at each joint, l!#iile 
resolving the member forces at any joint, in the released 
structure, the number of unknowns should be only three. 

Since the equilibrium equations are solved in a computer, 
the numbering of joints should be done in a particular 
sequence so as to meet the above requirement. For the 
present analysis of tower, by flexibility method, the joint 
numbering in the basket portion is done as shown in Fig. 4,1, 
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The numhering of joints In the tower “body, starting from 
40, is as shown in Mg, 2,7. 

The number of loading conditions to be considered 
for the analysis of the tower is 3 as concluded in Chapter 2. 

The time taken, on IBM 7044 computer, to analyse 
this tower for the three loading conditions using 
flexibility method is 37 seconds. 

The same tower, Pig, 2.7, has also been analysed 
by stiffness method. The number of degrees of freedom, 
n^, is 177. The half band width of stiffness matrix is 
30 for such double- circuit towers. This method takes 
only 31 seconds for the analysis, for the three loading 
conditions. The main core of the computer has been used 
for storing the stiffness matrix. 

It is also observed that stiffness method turns 
out to be better in computer memory requirements. However, - 
from the point of view of error accumulation, the flexibility 
method turns out to be superior. This is verified by 
comparing the resolved internal forces to exteinally 
applied loads. The comparison is presented in Table 4.1 
for typical joints. 

Since the stiffness method is superior both from 
the view point of computer memory as well as computer time, 
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Talol e 4.1 

Error Accumulated in 

Plexihilitv and 

Stiffness 



Methods 

of Analysis 




Joint TJo* 
(Pig. 2. 7) 


Plexihility Method 

Stiffness Method 


Is Pj 

UV 1 

iS 

Is Pj 

Is Pyl 

Iz P^' 

1 

0.00 

0.00 

0.00 

0.01 

0.01 

0.00 

6 

0.00 

0.00 

0.00 

0.01 

0.02 

0.01 

12 

0.00 

0.00 

0.00 

0.00 

0.05 

0.01 

20 

0.00 

0.00 

0.00 

0.02 

0.01 

0.00 

2 $ 

0.00 

0.00 

0.00 

0.00 

0.03 

0.00 

34 

0.00 

0.00 

0.00 

0.01 

0.01 

0.00 

40 

0.00 

0.00 

0.00 

0.00 

0.02 

0.01 

44 

0.00 

0.00 

0.00 

0.00 

0.01 

0,00 

48 

0.00 

0.00 

0.00 

0.00 

0.01 

0.00 

52 

0.00 

0.00 

0.00 

0,00 

0.00 

0.00 

56 

0.00 

0.00 

0.00 

0*00 

0.01 

o 

o 

-9 

o 


TalDle 4.2 Allowa'ble Maximujn Slenderness Ratio 


SI. No. 

Type of Memder 


Limiting L/'r 
Value 

1. 

Leg members 


150 

2. 

Members carrying computed stresses 

200 

3. 

Redundant members and those 
nominal stresses 

c arrying 

250 

4. 

Members subjected to axial 
only 

tension 

500 
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it is deemed to 'be more suitable for embedment in the 
automated optimum design of transmission line tower. 

Error accumulated in the stiffness method of analysis 
is negligible for all practical purposes, 

4- *5 Design of Members 

Sectional properties of available angle sections 

56 

as per Indian Standard Institution , commonly used for 
towers, are stored in the memory of the computer in 
ascending order of their areas. The first 54- are equal 
sections and the next 35 are unequal angles. The other 
sectional properties stored are, (1) least radius of 
gyration, (2) radius of gyration about x«-x axis^^ 

(parallel to smaller of the two sides) of the unequal 
angles, r^, (3) width and thickness of angles, and (4) 
width of connected side of the unequal angles. 

Members of the tower are put into a certain number 
of groups, as described in configuration generation. 
Section 2, 3.2.1, which is also useful for the design of 
members, Eor example all the leg members in a panel form 
one group and all the diagonal members in the same panel 
form another group. 

After the analysis is carried out, the member 
forces are compared to get the maximum and minimum forces 
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among eacTi group of members. These when multiplied by the 
appropriate factor of safety, presented in Section 2,4,2, 
gives the design forces of the group. The list of equal 
angles is then scanned (for diagonal braces, unequal angles 
starting from 55 in the list) from the beginning. The 
first section is chosen as the designed section of the 
group which satisfies the code specifications , They are 
listed below, 

(1) Compressive stress in the member should not exceed 
the permissible value, given in JSippendix 2, 

p 

(2) Tensile stress should not exceed 2600 kg/ cm , 

( 3) Slenderness ratio of the member adopted not to 
exceed the values recommended. Table 4.2, 

( 4) Thickness of angle selected for the members should 
not be less than 4 mm. 

The total weight of all the members, designed to 
meet the above requirements, is the weight of the primary 
members. The weight of all secondary braces which has been 
computed earlier. Section 2, 3. 2. 2, is then added to get 
the total weight of the tower, 

4.4 Results of Optimization 

The typical double- circuit tower, Pig, 1.1, subjected 
to three different critical loading conditions, discussed 
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earlier (Pig. 2.6), is considered for tlie optimization 
study, The following parameters have been kept constant 
throughout the study, 

(1) Plan of basket is 1.2m square and is constant 
throughout as shown in Pig. 2,6. The height of 
panels in the basket are not altered, 

(2) The base of the tower is square. 

(5) Tb.e weight of secondary braces in the basket is 

a constant, since the geometry does not change in 
this portion of the tower. Hence the weight of 
these secondary braces is not computed. Thus the 
weight of the tower indicated on tables and figures 
does not include the constant weight of secondary 
braces in the basket portion. 

The optimization study is carried out with different 
starting points to ascertain whether the local minimum is 
a global minimum. Studies are carried out on 15m, 20m, 
and 25m tower bodies. The results are presented in Table 4.3 
to 4,5 respectively. In these tables the panel heights 
are in metre and weight in kilogram. The average time 
taken for one function evaluation (weight of tower) varies 
from 11,0 to 15.6 seconds on IM 7044 computer. The CPU 
time for individual cases are given in the tables. 
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The reduction in the weight of the tower at the 
end of each one-dimensicnal minimization is plotted in, 

Pig. 4.2, for a few representative cases, A close look, 
at the results presented in Talle 4.3 to 4,5 and in Pig. 4, 2 
reveal the following: 

(1) The variation in the base width of tower has 

marked effect on the weight of the tower as compared to the 
variation in panel heights. This observation can be made 
from the graphs in Pig. 4.2, in vdiich iteration 1 and 
corresponds to base width. The reduction in the weight of 
the tower in these iterations, i.e., due to change in base 
width, is appreciable as compared to change in panel heights. 
This is explainable from the fact that increased base width 
leads to long members both for the leg as well as for the 
diagonals. In that case, even though the forces in these 
members are smaller, the increased length of these members 
lead to heavier sections. 

Moreover, for the diagonal members, it is observed 
that these are, in general ^ governed by the slenderness 
ratio restrictions. Therefore, the compromise between 
reduced force due to increasedbase width and reduced length 
of members due to reduced base width decides the base width 
corresponding to the minimum weight of the tower. Prom the 
optimization study, it is observed that, for the minimum 


Weight kg 
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weight tower, the base width turns out to be around 4 -ei for 
body heights of l5m to 25m. 

(2) The study further reveals (refer Tables) that for 

the minimum weight towers, number of panels should be 
chosen such that the height of the central panel(s) of the 
body is around 2,5m, Por an optimxhn configuration the 
panel heights are to be chosen such that the top most is 
the smallest and the bottom most has the largest panel 
height. Optimal panel height turns out to be around 2m 
and 5m respectively for top and bottom most panels. The 
intermediate panel heights may be varied linearly. This 
gets confirmed from the results of 25m body tower with 
10 panels given in Table 4.5 (columns 6 and 7) wbere 
initial design was taken close to the optimum. 

( 5) To study the changes in stiffness of the tower, 

the area of leg members at the starting and optimum designs 
are plotted in kig, 4.3. The plot correspond to the three 
towers for which optimum design results are presented in 
Fig. 4.2, The leg member areas are plotted against the 
mid-height of panels above ground. It can be observed 
from this figure that the changes in the member areas in 
the basket portion are negligible. However, the area of 
leg members in the body change considerably depending upon 
the height of tower. For the tower of 15m body height, the 




I 



j. 4*3 Area of leg members along the height of tower 1 
(Static loading) | 
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area of leg memlDers in the body increases while the base 
width decreases (columns 4- and 5, Table 4.3) at the 
optimum design. This is not so as the height of tower 
increases. Eor the tower of 25m body height, the area of 
leg members decreases while the base vddth remain almost 
the same (columns 4 and 5, Table 4.5). This shows that 
the minimum weight tower is tending to be more flexible 
as the height of the tower increases. Implicit in this 
approach of study is the fact that bending stiffness of 
the tower is primarily dependent upon the stiffness of its 
leg members. 

The above observation warrants the need to check 
the stability of high tower. It is only this study which 
shall indicate whether minimum weight obtained is at the 
cost of overall stability of the tower. Hence the stability 
of the minimum weight tower with 25 m body height has been 
checked using truss- element geometric stiffness matrix 
(initial stress matrix), Kq. The total stiffness matrix 
of the tower is obtained by superimposing Kq. and the 
conventional stiffness matrix, K. The details of the 
method are described in Chapter 15 of Ref. 83, After 
elimination of all lower triangular elements of the total 
stiffness matrix, the diagonal elements are found to be 
of the order of 10^, This shows that the determinant S K+K^l 
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Table 4-. 5 Results oi 25m Body Tower (Static Loading) 
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is very high (>> 1) and hence the tower is quite stable, 
Therefore, it is not necessary to embed stability consideration 
in the framework of optimum design of thjse towers. However 
a final check on the overall stability of the minimum 
weight tower is still recommended, 

(4) It is also observed, from the function values 
(weight of tower), that the reduction in weight after 
first cycle of minimization is negligible. The graph for 
15m body in Fig. 4,2 gives the function values through 'two 
cycles. A similar tendency is noted for higher towers also, 

So it is good enough to stop the optimization process after 
moving through the first pattern direction, 

(5) The results presented in Table 4.3 to 4.5 gives 
an idea of convergence to local minima, ALl the three 
designs, tabulated for 15ni body tower with 6 panels, havq 
converged to the minimum weight of around 3550 kg from three 
different starting points. Similar trend is noticed for 
20m and 25 m tower bodies, Hrom these it can be concluded 
with some degree of confidence that the minimum weight 
reported herein corresponds to the global minimum. 

^•5 Recommended Rational Sections 

From the angles picked up for the members of the 
tower, for all the results presented in the previous section, 



it is observed that certain angles are always avoided from 
selection, TMs raised the suspicion that some of the 
available angles are not rational. A look into the table 
of sectional properties reveal one of the possible reasons: 
when the angles are listed in the ascending order of their 
areas (Table A2.1), the other sectional properties, especially 
the least moment of inertia, do not follow the same 
distribution as that of the area. Based on this observation 
a study has been conducted. The study resulted in the 
reccmimendation of rational angles for the members of any 
truss structure. The detailed discussion and the list of 
recommended angle sections are presented in Appendix 2, 

Bor the recommended rational angles, the area as 
well as the least moment of inertia have stepwise increase 
in their values. Using these recommended angles, in place 
of the available I.S. angle sections, the optimum design 
of transmission line towers is carried out. The starting 
configuration is chosen based on the results obtained 
using available angles, in the previous section. The results 
of this study are presented in Table 4.6 for 15m, 20m, and 
25m body of the tower, A further reduction of 2,J5 percent 
is achieved for 20m and 25 m body heights by using suggested 
rational sections. However the reduction in weight for 
15m tower body is negligible, i.e.. 


minimum weights obtained 
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using I.S, angles as well as the recommended angles are 
nearly the same. This is due to the fact that for small 
area range, the rational angles recommended are also available 
in ISI Handbook (refer to areas in Table A2.1 and A2.2), 
which are required for 15m tower body. This is not so for 
large area range and hence the reduction in weight of the 
tower is considerable for 20m and 25m body heights. 

Results of Table 4.6 further reveal that the . 
optimum base widths of tower obtained are less than those 
obtained with available angles. The optimum base widths 
obtained vary from 3.60m to 3.90m for the range of heights 
considered in the present work. There is no change in the 
conclusions arrived already in respect of panel heights. 



CHAPTER 5 


OPTIMUM DESIGN POR DYNAMIC LOADS 
5 .1 Introdu ction 

Structural design for d 3 mamic loads consists of 
three steps, viz,, (1) idealization of structural loading 
as a function of time, (2) study of structural characteristics 
via modal analysis, and (5) response prediction inorder to 
evolve a safe design. The first step of load idealization 
is already described in Chapter 2 , The present chapter 
starts with modal analysis of transmission line tower 
followed by the description of response prediction through 
modal superposition method, Optimimi configuration design 
under dynamic loads is then carried out. The results of 
dynamic analysis and optimum design are tabulated and 
discussed in this chapter, 

5 *2 Modal Analysis 

The eq^uation of motion for a system of structural 
elements under dynamic loading can be written as 

Mu + Cu+Ku = p (5.1) 

where M, C and K are respectively mass, damping and stiffness 

^ -J ^ 

matrices of the system; u , u » n and p are the displacement 
velocity, acceleration and load vectors respectively. In 



96 


structural application, the K and M matrices are the result 
of the finite element idealization®^ of the actual structure. 
The size of matrices involved is n^ x n^ where n^ is the 
number of degrees of freedom of the system. 

The equation of motion for the free vibration of 
the structure is obtained from Bq. (5»1) by setting the 
right hand side as zero vector and assuming the system 
undamped as 

M u + K u = 0 (5 .2) 

Substituting for the displacements, 

u = 0 sinw (t - tg) , (5.3) 

we obtain the generalized eigenvalue problem 

K 0 = 0)^ M ? (5.4) 

The eigenvalues give the natural frequencies 
of the system and the eigenvectors are the corresponding 
modes of vibration. The complete solution to Eq. (5.4) 
can be written as 

K S = Q® M $ (5.5) 

in which the columns in ffi are the eigenvectors an^ 

= diag,( ( 0 ?) . 

The basis is now changed from the physical 
coordinate system which was used to establish Bq, (5.1) to 
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the M-orthonormal "basis of eigenvectors. Using 

u = ® q (5.6) 

where the vector, q , list the coordinates in the new basis, 
The equilibrium equations, Eq. (5.1), become 

q + Al + q = p (5.7) 

T 

in which A = i 0 i and is assumed to be diagonal « This 

86 

requires that the damping matrix is of a restricted form , 
Eq. (5.7) consists of n^ uncoupled equations which can 

QT^QQ 

readily be solved ^ , 

5.2.1 Eigenvalue Problem 

The most time consuming step in the dynamic analysis 
is the solution of the eigenvalue problem, Eq, (5 .5) . If 
the order of the matrices is large, the computer time 
required to solve for all the eigenvalues and vectors 
can be enormous. However, the structure may respond 
primarily in only a few modes from the lowest end of the 
spectrum and the contribution of the other modes may be 
negligible. In this case a method of solution is needed 
which calculates only the required frequencies and mode 
shapes economically. 

The detailed theory on matrix eigenvalue problem 
is available in Ref. 90. However, in structural applications 
the approach is slightly different. Based on the structural 
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stiffness properties, minimum energy principle, and from 
the experience gained to date, two approaches emerge for 
partial eigensolution of large systems, viz., 

(1) Pol 5 momial iteration method 

(2) Suhspace iteration method. 

In the first method, the zeros of the characteristic 
polynomial, p(X), given by 

p(X) = Ik - XM| (5,8) 

is searched in an iterative manner. Hence this method is 

advantageous only for systems having small hand width^'*’”^^ « 

On the other hand, the suhspace iteration method, can he 

94 95 

used effectively for large handed systems ’ , This 

method uses the Rayleigh minimum principle. Most of the 
methods for partial eigensolution can he shown to fall in 
the latter category. 

for the size of matrices involved in the present 
study, discussed in Chapter 2, suhspace iteration method 
is considered most suitable , The details are described in 
the following section. 

5 .2 ,2 Meth od of Solution of Eigenvalue Problem: 

Suhspace Iteration Method 

96 

The suhspace iteration method , in essence, consists 
of the following three steps: 
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(1) Establish, q starting iteration vectors, q >; p, where p 
is the number of eigenvalues and vectors to be 
obtained . 

(2) Use simultaneous inverse iteration on the q vectors 
and Ritz analysis to extract the ’best’ eigenvalue 
and eigenvector approximations from the q iteration 
vectors , 

(3) After convergence to the subspace, use the Sturm 
sequence check to verify that the required eigenvalues 
and corresponding eigenvectors have been extracted. 

The solution procedure has been called subspace ' 

iteration method, because the iteration is equivalent to i 

iterating with a q- dimensional subspace and should not be | 

r 

regarded as a simultaneous iteration with q individual S 

vectors. [ 

h 

I 

The basic objective in the subspace iteration i 

t 

method is to solve for the lowest p eigenvalues and the ! 

i 

corresponding eigenvectors satisfying [ 

I 

K S = A M s (5.9) I 

where A = diag( X^) and = 10.^,....., 0p3 • 

i 

The essential idea of the subspace iteration 
method uses the fact that the eigenvectors in Eq» (5.9) 
form an M- orthonormal basis of the p- dimensional least- 
dominant subspace of the operators K and M, which is called 
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. In the solution, the iteration with p linearly 
independent vectors can therefore he regarded as an iteration 
with a suhspace* Ihe starting iteration vectors span 
and iteration continues until the converged solution 
theoretically spans The fact that iteration is performed 

with a subspace has some important consequences* The total 
number of required iterations depends on how 'close* Ej is 
to and not on how close each iteration vector is to an 
eigenvector* Hence the effectiveness of the algorithm lies 
in that, it is much easier to establish a p- dimensional 
starting subspace which is close to than to find p vectors 
where each is close to a corresponding eigenvector* Also, 
because iteration is performed with a subspace, convergence 
of the subspace is all that is required and not of individual 
iteration vectors to eigenvectors* In other words, if the 
iteration vectors are linear combinations of the required 
eigenvectors, the solution algorithm converges in one step* 

The algorithm presented below is the procedure 
of step(2)of subspace iteration method. 

For k = 1,2, *.**, iterate from to Ej^^^ as follows: 

(1) Carry out the inverse iteration, 

K = M (5.10) 

where X| is the matrix consisting of the set of 
assumed starting iteration vectors and is the 
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intermediate improved set of vectors for Ritz 
analysts. The constraint on the choice of is 
that none of its vectors he orthogonal to the 
corresponding eigenvectors , 

(2) Find the projections of the operators K and M onto 

®k+i = 


^ta-1 

II 


(5.11) 

^TtC+l 

= Vi 

M X^_^^ 

(5.12) 


(3) Find the eigensolution of the projected operations, 
i.e., solve 



^k+1 

Vi Vi Vi 

(5.13) 

Find 

an improved approximation to 

the eigenvectors 



V-1 Vi 

(5.14) 

Now 

\+i 

^ ^ as h 


and 


— S } 



In the suhspace iteration, it is implied that the iteration 
vectors are ordered in an appropriate way, i.e,, the iteration 
vectors 0^ , 02 » ••••» ar© stored as the first, second,,,.., 
columns of . 

The effectiveness of the algorithm lies to a large 
extent in that a few iterations can already give good 
approximations to the required eigenpair. The reason is 
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that one suhspace iteration in Eq., (5.10) to Eq, (5.14) is> 
in fact, a Pdtz analysis described in Appendix 3. In other 
words, the eigenvalues from lowest end of the spectrum are 
approximated best in the iteration and all eigenvalue 
approximations are upper bounds^"^’^® on the actual eigenvalues 
sought . 

A notable observation in the subspace iteration 
method described above is that and in Eq, (5.1l) 

and Eq. (5.12) respectively, tend towards the diagonal form 
as the number of iterations increases; i.e., and 

are diagonal matrices when the columns in store multiples 

of eigenvectors. Hence the generalized Jacobi method, 
presented in Appendix 4» is used very effectively for the 
solution of the eigenproblem in Eq. (5.15) . Ihe size of 
matrices in Eq. (5.15) is only p x p , 

An important aspect is the convergence of the method. 
Assuming that in the iteration the vectors in are ordered 

in such a way that the i-th diagonal element in is 

larger than the (i-1)th element, i=2, ...»p, the i-th column 
in converges linearly to and the convergence rate is 

. Although this is an asymptotic convergence rate, 
it indicates that the smallest eigenvalues converge fastest. 
In addition, a -higher convergence rate can be obtained by 
using q iteration vectors, with q ^ p. However, using more 
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itieration vecijors will only increase the computer effori 
for one iteration. 

It is apparent that the closeness of to ^ is 
the measure of convergence which can he quantified hy the 
computed eigenvalues (or eigenvectors), let the computed 
eigenvalues in the iterations (k>-1) and k respectively he 
and ^ then the ratio 

is used as a measure of convergence. For example, for the 
eigenvalues to he accurate to about 4 digits, it is necessary 
to iterate until this 8 is less than 10 . 

The resulting eigenvalues are checked hy using 

the Stum sequence property^^ . This property states that 

T 

in Gauss elimination to evaluate m = K - dM, in which 

y is the shift and 1 is the lower triangular matrix, the 

number of negative elements in the diagonal matrix 1 is the 

number of eigenvalues smaller than U , In order to use a 

meaningful v , it is necessary to find hounds for the exact 

eigenvalues using the calculated values • A 

conservative estimate for a region in which the exact 

94 

eigenvalues lie is given hy 

0,994^^ < \ < 1.01 (5.16) 

Further details of the method qnd the computational time 
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estimates are available in ^ef, 96. 

5 .2 .3 Eigen solut'' on of Transmiss i on Line lower 
5 .2,3.1 Starting Sutsp ace 

The first step in the subspace iteration method 
is the selection of the starting iteration vectors, . 

The starting vectors for the present problem of transmission 
line tower is obtained from static deflection shapes. 

However, for the general structural application, the 
selection of starting vectors can be chosen based on the 
elements of K and M matrices « 

The transmission line tower has been modelled as 
n^ degrees of freedom system (n^ = 3n. where n is number of 
joints, above the fixed base of the tower) . Thus a complete 
eigensolution will consists of eigenvalues and associated 
mode shapes . However, a reasonable accurate d 5 mamic analysis 
of the system can be performed by estimating only first few 
eigenvalues and associated eigenvectors from the lowest 
end of the spectium. Por this purpose the tower can be 
thought of as a tapered cantilever beam in space having 
almost equal stiffness in the transverse and longitudinal 
direction. The actual variation in stiffness along these 
two direction is due to cross arms. Such a tapered cantilever 
beam is capable of having principal modes of vibration along 
the transverse direction, longitudinal direction and about 
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c entroidal axis. The principal modes of vibration along 
transverse and longitudinal directions are bending inodes 
while those about the centroidal axis the torsional modes. 
The first two principal modes along the transverse direction 
as well as about the centroidal axis for a cantilever beam 
are shown in Pig, 5.1. The corresponding static deflection 
curves for the transmission line towers are obtained by 
applying suitable static loads as follows: 

The first mode of Pig. 5.1(b) is approximated 

by applying a concentrated load at the tip of the tower 

along the transverse direction. The first mode in the 

longitudinal direction is approximated, in a similar way, 

by applying a concentrated load at the tip of the tower in 

the longitudinal direction. The second mode of Pig, 5. 1(b) 

is once again approximated by applying two concentrated 

loads one at the tip of the tower and the other 3 times in 
99 

magnitude of that at the tip, applied at the level of 
bottom cross arm of the tower, in the opposite direction. 

The first mode of Pig, 5.1(c) is approximated by applying 
a couple at the top conductor points of the tower. The 
second mode of Pig. 5.1(c) is approximated by applying two 
couples, of same magnitude, one at the top conductor points 
and the other at the bottom conductor points in the opposite 
direction. All these static loads are shown in Pig, 5.2, 
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(b) Bending modes in x-y plane 



(c) Torsional modes 


Fig. 5-1 Typical mode shapes of a cantilever beam ' 




Fig52 Loading conditions j 

[ Defieqtcd shape corresponding to toad 
condition as starting vector to 
extract j-th mode shape] . j 

f 

i 
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The starting iteration vectors, obtained from 
static deflection shapes, are to be arranged in such a way 
to get the eigenpair convergence in the increasing order 
of their eigenvalues , This can be easily achieved by trial 
and error method or by solving a series of individual 
problems and arranging the eigenvalue solutions in ascending 
sequence* 

Six eigenmodes in the increasing order of eigenvalues 
as listed below have been extr^acted in the present work. 

1 , First bending mode, in the transverse direction . 

2 , First bending mode in the longitudinal direction 

3 , First torsional mode 

4 , Second torsional mode 

5 ^ Second bending mode in the transverse direction 

6, Second bending mdde in the longitudinal direction. 

5, 2. 3. 2 Results 

A tower of 15m body height is chosen for the study 
of eigenpair. The base width of tower is 5»4m. Number of 
panels in the body has been taken as 6 , The size of the 
stiffness matrix is 177 square and its half band width is 30. 
Lumped mass formulation is used and hence mass matrix turns 
out to be diagonal. 

The masses lumped at the joints are very small and 
this posed computational problem. To overcome this problem 
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scaling has been applied for the elements of the mass matrix 
as described below: 

Consider the eigenvalue problem 

K 0 = XM 0 (5 .17) 

Rewriting the above equation, with the introduction of a 
constant c, we get, 

K ? = (0 M) 0 

c 

i.e., K ^ = "Xm 0 (5.18) 

The eigensolution of Eq. (5.18) can be obtained without 
any computational difficulty, since the elements of M are 
c M , c being any suitable constant. The eigenvectors, that 
are obtained from the solution <of Eq. (5.18) , do not change 
where as the eigenvalues so obtained have to be multiplied 
by c to get the eigenvalues of Eq. (5.17) . 

To get a suitable value for the constant c, the 
eigenproblem of the transmission tower, under discussion, 
has been studied with different values for c. The results 
of this study are presented in Table 5.1(a) , (b) and (c) . 

Erom these tables, it is clear that the scaling 
constant of 1000 and 10000 have pruduced equally good 
convergence, whereas it is slow with c = 100 especially for 
higher eigenvalues . So the scaling constant of 1000 has 
been used for the eigenpair extraction. 
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T abl e 5 . 1 ( a) Resu lts of Eigen^- Iteration with Scaling Paotor 
c = 100 

V. -m r iftwi m 


Iteration ipproximation to Eigenvalues of the given problem, 
NO. Eq*(5*l7) 



1 

2 

3 

4 ' 

5 

6 

1 

844*68 

874.08 

1619 .36 

7319.56 

7962 .39 

9063.24 

2 

844.66 

874.06 

1617.43 

7125 .40 

7957 .30 

9051 .26 

3 

844.65 

874.06 

1617.28 

7104.52 

7957.20 

9051 .02 

4 

844 .66 

874.06 

1617.42 

7104.19 

7957.20 

9051 .13 

5 

844.66 

874.06 

1617.41 

7103.90 

7957.17 

9051 .13 

6 

844.66 

874.06 

1617.33 

7103.84 

79 57.96 

9051 .13 

15 

844.66 

874.06 

1617.42 

7103.88 

7957.10 

9051 .19 

16 

844.66 

874 .06 

1617.41 

7103.79 

7957.07 

9051 .23 


Table 5 * 1 ( b) Resul ts of Eigen- Iteration with Scaling Factor c=1000 


Iteration Approximation to Eigenvalues of the given problem, 

No. Eq. (5.17) 


1 

2 

3 

4 

5 

6 

844.69 

874.07 

1619.36 

7319.56 

7962 .44 

9063.14 

844.66 

874.06 

1617.43 

7125.46 

7957.30 

9051 .26 

844 .66 

874 .06 

1617.42 

7106.34 

7957.21 

9051 .13 

844.66 

874.06 

1617.42 

7104.19 

79 57.21 

9051 .13 

844 .66 

874.06 

1617.42 

7103.90 

7957.21 

9051 .13 

844 .66 

874.06 

I 6 I 7.42 

7103.86 

7957.21 

9051 .13 


6 
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Table 5.1(c) Resul ts of Eigen-Iteration with Scaling factor 
c = 10000 


Iteration Approximation to Eigenvalues of given problem, 

No. Eq. ( 5 . 17 ) 



1 

2 

3 

4 

5 

6 

1 

844.69 

874 ,08 

1619.36 

7319.56 

7962.87 

9061.61 

2 

844.67 

874 .06 

1617.44 

7125.46 

7957.30 

9051 .26 

3 

844.66 

874.06 

I 6 I 7.42 

7106.34 

7957.21 

9051 .13 

4 

844 .66 

874 .06 

I 6 I 7.42 

7104.19 

7957.20 

9051 .13 

5 

844.66 

874 .06 

161 7 . 42 

7103.90 

7957.21 

9051.13 

6 

844.66 

874.06 

1617.42 

7103.86 

7957.20 

9051 .13 
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Using the subspace iteration method, eigenpairs 
of a transmission line tov^er are computed. The tower 
considered has 15m body height and 4.33m base width. Ther^ 
are 6 panels in the body. The results of eigenvalue 
iteration are presented in Table 5.2. 

Table 5 .2 Eigenvalues of Transmission Line Tower 


Iteration 

No . — 


Eigenvalues 




1 

2 

3 

4 

5 

6 

1 

691 .39 

710,83 

1230.85 

6171 .59 

7046.35 

8041 .77 

2 

691.39 

710.82 

1230.15 

6014'.02 

7042 ,86 

8033.62 

3 

691 .39 

710.83 

1230.15 

6002 .25 

7042'.82 

8033.55 

4 

691 .39 

710.83 

1230.15 

6001 .26 

7042 .81 

8033.55 

5 

691 .39 

710.83 

1230.15 

6001 .16 

7042 .81 

80 33.55 

6 

691 .39 

710.83 

1230.15 

6001 .15 

7042 .81 

8033 .55 


The mode shapes corresponding to the first six 
eigenvalues listed in Table 5.2 are sketched in Mgs, 5.3(a) 
(T) . In these figures, the bending modes have components 
only in a plane, i.e., either in X-T plane or Y-Z plane. The 
components in the perpendicular plane is negligibly small 
(of the order of lO"^ or less Pereas those in the plane is 
of the order of 10 ) . ' This confirms that the eigenspectrum 

is a series of planar bending modes and torsional modes 
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Starting from the lowest end of the spectrum, 

The lowest eigenvalue is 69 1 ,39 and hence the 
fundamental natural frequency m^=Y691,59 = 26.3 rad/sec. 

The large values of are due to the fact that the elements 
of mass matrix are very small in comparison to the elements 
of the stiffness matrix. 

The second eigenvalue is very close to that of 
the first. Theae two eigenvalues correspond to the first 
tending modes in the mutually perpendicular planes, Pig, 5. 3(a) 
and (h) . The stiffness of the tower is almost the same 
in these two planes. The small difference is only due to 
the presence of cross arms in one of them and hence, the 
first two eigenvalues turns out to he close . 

The two torsional modes lie adjacent to one another 
in the series. The eigenvalues of these modes are very 
well spaced. Prom the eigenvalues given in Table 5 '.2, it 
is observed that the first three eigenvalues, corresponding 
to the first mode in each of the principal directions, being 
close to each other can be called a set. The rest three 
eigenvalues belonging to the second set of eigenmodes are 
very well separated from the first three, 

5 .2 ,4 Parametric Stud y 

A study on the variation of eigenvalues with height 
of tower has been conducted. Eigenvalues are extracted for 
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transmission line towers having total height of 28.0m, 30,5m, 
33 «0m, 35.5m and 38,0m, The height of basket is taken as 
13m (Pig, 2,6) and is constant for all the towers considered 
in this study. The width of base has been fixed at 4.33m 
for all towers based on the optimum design results for 
static loads. The panel heights are also decided from 
the observations on the optimum design results, Chapter 4, 

The results of parametric study are presented in Table 5.3. 

The base width being kept constant, tiie results 

presented in Table 5.3 are essentially a function of tower 

height. These results are plotted in Pig. 5.4(a) and 

Pig. 5.4(b) against the height of tower. The results are 

2 

drawn in two separate figures, since is about 10 times 

2 

of co-j , i.e,, spaced much apart. 

As discussed in the previous section, the first 

tvro eigenvalues are close to each other and the difference 

goes on diminishing for increasing height. Beyond certain 

2 

height, say 40m (Pig, 5, 4(a)), the eigenvalues and 

2 

Wg of the space transmission tower become equal, but the 

eigenmodes are different. This shows that the effect of 

cross arm on overall stiffness of the tower becomes negligible 

2 

for high towers. The third eigenvalue , corresponding to 
torsional mode, is separated apart and this does not come 
close to other lower bending modes for all heights of 
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Fig.5*4(b) Height of tower versus ouj (i=4,5,6) 


125 


transmission line towers considered, 

2 

In Pig, 5.4(13), 0 )^ correspond to the second torsional 

2 2 
mode, to the second lending mode in X-Y plane and Wg 

to second bending mode in Y-Z plane. This order of 

eigenmodes remained the same (also listed in page 108) for 

heights upto 57m, The order gets changed for higher 

towers, i,e,, for tower height greater than 57m, Pig. 5.4(b). 

In the changed spectrum, the fourth and fifth eigenmodes 

2 

interchange their position. The absolute slope of to^ curve , 

ft 

Pig. 5, 4(a), which corresponds to the first torsional mode, 

2 2 

fast increases than that of and», 0)2 curves as the 

height increases , Furthermore, the absolute value of 

2 2 2 
slope of the curve is more than that of and for 

the range of heights considered. These observations 

show that the eigenvalues corresponding to bending modes 

approach zero faster than that of torsional modes a,s the 

height increases. Hence, for high towers, the torsional 

modes move away leaving the bending modes to dominate the 

lowest end of the eigenspectrum. 

For tall towers, the torsional modes are driven 
away from the lov/est end and hence for slender towers only 
bending modes need be extracted for a reasonably accurate 
dynamic analysis. Furthermore, as noted earlier, for 
slender towers, the eigenvalues corresponding to the 
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"bending modes in two mutually perpendicular* planes become 
equal. Thsse two observations state that the modelling 
of structure can be further simplified and idealized as a 
plane frame. This boils down to the conventional d 5 mamic 
analysis of multistoreyed buildings and television towers, 
wherein the adjacent mass points are idealized to lie in 
a vertical axis and are assumed connected in between by a 
single equivalent plane frame member. 


5,3 Dynamic Response 

5,3.1 liyin d Load Analysis 

The n^ uncoupled equations, Eq, (5.7), can be 
solved if the diagonal matrix, A , derived from damping 
matrix C is known. But in practical dynamic analysis the 
matrix € is not known. However it is usual practice to 
adopt an equivalent fraction of critical damping, for the 
structure considered. The same has been adopted in the 
present work. Therefore, the j-th equation of the uncoupled 
system, Bq, ( 5 . 7)5 l.e .5 j-th modal equation of motion, 
turns out as 

2 

*qyt) + 2p . 4 .(t) + 0 ) . q.(t) 0^^ p^(t) (5.19) 

where suffix i denotes i-th degree of freedom, 

p . = K is the fraction of critical damping 

i-th component of the j-th mode shape vector, which 
^ 3 


are M- orthonormal 



Pj^(t) - load function at i^th degree of freedom. 


Assuming initial conditions to "be zero, the solution 
of Eq_, (5.19) is given®® as 

^ t 


qj(t) ^ m^^(t^T) dx} 

( 5 . 20 ) 


where m 


D<i 


_ fa -2 2 / _2 ■ 

- h- . 


The functions p^(t) can he replaced, for computational 


purposes, as a product, 

Pj_(t) = fj_(t) 


( 5 . 21 ) 


where is a constant, i.e,, maximum value of load Pj_(t) 
over the time interval considered, f^(t) is a function of 
time which takes a maximum value of + 1 .0 . Substituting 
Eq. (5,21) in Eq. (5.20), the response equations turn out 
to he, 




rif 

= S 
i=1 


^i: 


Pi{ 


0) 




t 

f fi('r) 


^-P3(t-T) 


sin 


W;j^(t- T) 


dx } 


( 5 . 22 ) 


Introducing the concept of dynamic load factor (DLP) , which 
is defined as the ratio of dynamic displacement to static 
displacement, the maximum value of q^—— ihe modal 
participation factor, is given hy 


I 


^d,max =.5^ 



(5.23) 


where ij defined as the maximum DIiF for the 

unit load fj_(t) due to oscillation under natural frequency 
rad/sec . 

The peak values of are being known from Eq.(5.23), 
the upper hound on displacement can he computed by 
superimposing the modal contributions . 

Uf 

%,max Ijjmax I (5.24) 


The summations in Eq. (5.20) through Eq. (5.24) shall go 
from 1 to s, only, when s out of n^ eigenpairs from the 
lowest end of the spectrum are used to get a reasonable 
approximate analysis. Hence, in the rest of this chapter, 
superposition of only first s modes are considered for 
presentation. 

In order to compute the maximum, stress at any point 
on the structure, we need to know its displacement mode. 

We have available to us two choices: (1) to use the upper 
bound displacement vector whose components are given by 
Eq. (5.24), or (2) to use the natural modes of the structure, 
0’^ , d = 1,2,,..,s and then superpose the resulting stresses 
corresponding to each natural mode in proportion to the peak 
modal participation factors, Eq.(5.23). 



The use of the upper bound displacement vector 
seems unreasonable because its components represent an 
upper bound to the corresponding peak displacement and 
is merely an approximation (though reasonable) to the 
actual displacement of the structure. Moreover, its 
use may lead to a conservative estimate of stress at one 
point on the structure and simultaneously may result into 
a nonconservative estimate of the stress at some other 
point on the structure. 

The use of natural modes as indicated in choice (2) 
seems more logical, since the structure in any case vibrates 
in its natural modes. We have, therefore, the upper bound 
stress on any member i as 

^i,max “ 1 ®ij ^j,ma;xl (5.25) 

where is the stress in member i due to mode shape 3 
taken as deflection vector, 

5, 3.1.1 D 3 mamic Load Factor (DLF) 

The wind loads on joints at peak level have been 
calculated and shown in Fig. 2,6, The time function, f (t) , 
is shown in Fig. 5.5. The method of computing (^hP)jjjg_ 2 c 
this unit load-time function is described below. 

Dynamic load factor (DIF) is given by, from Eq.,(5.22), 
for any general load function, f(t) as: 



DLP = 


03 

"oT 


/ f('^) ginw,(t-T) dT^ ( 5 . 26 ) 

d o 

The above e:Q3ression is nothing but the displacement 
function, y(t) , sines static maocimum displacement, yg-j-> 
of f(t) is 1 .0 « Hence, without any loss of generality, 
we Can write the above equation as 
0)^ t 

y(t) = ~~ / di; (5.27) 

The deteiministic wind load function, f(t), shown 
in Fig. 5.5 has four distinct segments as follows: 


(1) 

f(t) 

II 

= b = H 

0 < 

(2) 

f(t) 

= P 

0 

> 

S ^ S. "^2 

(3) 

f(t) 

= p+(1 . 

. t 

; t2 < t _< t 

(4) 

f(t) 

= 1 .0>-(1 

•o-p)t^ ;%=Vt 



(5 . 28 ) 


The rest of the load function is the cyclic repetition 
of load segments (2), (5) and (4) listed above. Substituting 
the above functions in Eq, (5.27), we get y(t) for different 
regions as given below: 

Region 1 : o ^ ‘t-j 




2 t 


/ 2— sin w,(t-T) dm 


‘^d ° ■*^d 

P 2p 

t t >- ' { 2P COS (■ 


0) 




10 ^) sin u^jt}] 
(5.29) 


(5.30) 


[1 




U. 


Region 2 : 't-] .< t _< t2 

y = sin cos o.^- 

Slna..( 

do ^ 


U),t) 


d-r 


p + e' 




:^o + P^O 


pp 


0 ) 

d 


y = [ Yq cos w^t 


2 

+ w (yo~p) 1 

^ sinw^tJ 


sin w^t ] 

(5.51) 

(5.32) 


from Eq^, (5.29) 
from Eq, 


y(t=t^) 

. (5.30) 


where y^ 

(t=t^) 

Region 3: t2 5, increasing from p to 1.0 

,VPyo . 

SIJI 10 ^-t + y^ COS lU^t) 

d 

2 ^ 

+ f [p+(1.0-p) •-™-] e' ^ ^ sin w^( 

d o tjj^ ^ 

-0t ,VPyo . 

(__„„ a)^t + y^ 


y = e 


-pt 


= e” 


cos “^t) 



(cos 

w ,t + ■ 

CL I 

1-p 

') [t- 

2p 

e-Pt 

+ ( — 
^d 

2 

CO 

oF 

1-P 

(t ^ 

2P 



00 

) + e I 


“ — sin 0) ^t) J 
0 Q- 




W j) si 


^d 

(A cos w^t + 3 sin 
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where A = 


B = 


( 1 -p) 2 p 

y - p + — V— 7-' 


(pf;; ‘"d) 

“d 0)2 0 )^ 


y©^ y(’fc=*t2’"'*^i) (5.3i) 

^0== y('fc=V‘^l) ]3(1.(5.32) 

1"-P w-Bt 

y = + ® [ (B w^^PA) cos w^t-(pB + Aw ^) sin w^t ] 

(5.54) 

4 • "^3 Ji ^ "^4 > load decreasing from 1 ,0 to p 

-6t ,^o + Pyo . + 4.N 

g ^ sin W(jt + y^ cos w^t) 


y = e 


CO 


-t- 


(0 


= 1.0 


d o 
(1-p) 


•{; 

/ [1.0- (1.0-p) -~p- ] sin w , (t'-T) dT 

n 1 1 


■(t 


2 P 


oF 


') + e"'^'*' (A cos CO ^t 4 - B sin w^^t) 

, (5.35) 

1 ‘-P, 2 P 

“2 


where A = y^, >- I- ( 7 —-) 


^■td 


CO 


B 


_ + ^yo pP ('■"p) (^^ ^ Wfl ) 

— 1.1 '■> -a.ti I— I .^r 'a 1^ > 


0 ) 


^d ‘^d 


yo" y (■t = tj*-t 2 ) from Eq. (5.33) 

yo= y (t = t^~t 2 ) from Eq, (5.34) 
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Fig.&6 DCF curves for undamped case 



(c) ^*> = 25 rad /sec 


Fig.5-6 (Continued) 



2% damping 
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5 .3.1 .2 Maximum DLF 


The curves in Figs. 5.6 and 5.7 show that the 
maximum value of DLF occurs during the period of gust'. 

This suggest that the (TiIiF)jj^g^ can he obtained by searching 
only that time interval during which gust appears. Since 
the maxim'cim of y (or DLF) correspond to f <0, Eq, (5.39), 
the length of time for which > 0 can be ignored from 
such. Thus computations to obtain (Dli^P)jj^ax considerably 
reduced. To get an overall variation of (DIP) as a 
function of u , (DIE)j^g^ is computed for w ranging from 
0 .2 to 4-0 rad/ sec at intervals of 0 .2 rad/sec and the 
values are plotted in Pig. 5,8. 


The respfijnse (or DIP) during gust is very much 
dependent on the initial conditions at the start of the 
gust. Purthemore, the function f(t) is multimodal. Hence, 
the (DltP)jjiax c^i^ves in Pig. 5.8 are not monotic . 

5 .3.2 Earthquake load Analysis 


The j-th modal equation of motion of the system 
subject to earthquake load is given by 


Jil(t) + 2 ^. 4^(t) + a)| g^(t) = . 


rgCt) 

i=1 



mi 


where ^gCt) is the prescribed suppoit acceleration and 
9.(^) is the relative modal displacement with 


respect to the 



p=o-o 



Fig. 5-8 (DLFl versus 
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support. Proceeding similar to that of wind analysis, the 
maximum value of q^. can be obtained, for the present 
analysis, as 



where is the acceleration spectrum value corresponding 

to the natural frequency rad/sec. Prom Eq, (2,2) we have 

Sa,3 = “h 8= P I ^c, Sa <5-^2) 

The quantities in the above equation have been defined 
in Chapter 2 and the value of is read from Pig. A1 ,1 
for appropriate natural period and damping of the structure. 
The rest of the precedure for calculation of response of 
the structure follows exactly as in the case of wind load 
analysis . 

5 .4 Results of Dynamic Analysis 
5*4.1 Wind load Analysis Results 

Transmission line towers having body heights of 
15in, 20m, and 25m and a constant base width of 4*2m are 
analysed for time dependent wind load. Static loads 
simultaneously acting on the tower are treated as a special 
case of dynamic loads, i,e,, = 1,0 for all 

The maximum dynamic stresses in the members of tower are 
computed by modal superposition method imder three critical 
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loading conditions discussed in Chapter 2, The results 
are presented in Tables 5.4- to 5.6. 

The first colujan in these tables gives the 
connectivity of the leg members of the tower from the 
top. The first seven members correspond to the leg members 
of the basket and the rest to the body of the tower. 

These have been taken as representative member forces for 
the purposes of present discussion'. The second coliunn 
gives the estimate of maximum d 5 mamic force in the 
corresponding member by considering the superposition 
of only first three modes of vibration from the lowest 
end of the spectium in d 3 mamic analysis. The third column 
shows similar results considering six modes of vibration. 

The last column gives the corresponding forces due to 
static analysis of the tower, discussed in Chapter 4, 
considering maximum wind loads as static loads. The 
member forces in colunms 2, 3 and 4 are suitably modified 
by the factor of safety discussed in Chapter 2 to get the 
design forces , 

Comparison of values in columns 2 and 3 reveals 
that the member forces vary substantially in the leg members 
at the top while those at the bottom of the tower differ 
by about 10 percent'. The dynamic forces in the leg members 
of the body of the tower seen to be converging faster than 
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those in the basket'. Since the intent of the present work 
is to study the optimum configuration of the tower subject 
to dynamic wind loads and since the configuration of the 
basket portion is fixed from non- structural considerations 
it is reasonable from the point of view of accuracy as 
well as computer time to carry out the dynamic analysis 
by modal superposition of only first six modes from the 
lowest end of the spectrum. Consideration of hi^er modes 
shall increase the member forces only marginally. 

Once again from the comparison of results in 
columns 3 and 4 of these tables, it is observed that the 
member forces due to dynamic considerations are 9 to 
33 percent more than the corresponding forces due to 
static considerations. The percentage increase in member 
forces is more for higher towers. This confirms the earlier 
observation that the equivalent static wind load analysis 
is more conservative for small towers as compared to high 
towers. Therefore, to get a realistic assessment of member 
forces in transmission line towers a rigorous dynamic analysis 
is warranted. 
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TalDle 5 *4 M emlper forces in 15ni Body Tower ( 6 Pamiels) 


Joints of one of Dynamic analysis 

the leg members 3 jj^ode ““ 'TmodT””'” 

(Fig ^°2 superposition superposition 

Static 

analysis 

2 , 

8 

2.184 

3.214 

7.220 

8, 

12 

5.648 

7.194 

10 .166 

12, 

16 

13.617 

14.995 

15.904 

16, 

22 

17.313 

19 .106 

20 .855 

22, 

26 

21 .854 

27 .209 

26.241 

26, 

30 

28.290 

35 .671 

34.325 

30, 

36 

32.659 

38 .850 

36.607 

36, 

40 

40 .405 

44.764 

43.555 

40 , 

44 

45.336 

46 .257 

49 .086 

44, 

48 

47.034 

48 .375 

50 .289 

00 

52 

50.407 

53 ' .037 

53.215 

52 , 

56 

51 .929 

55.841 

53.939 

56, 

60 

54.189 

59.132 

55.710 


Forces in metric ton 
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Table 5 *5 Mjfflbe r Eorces in 2Qni Body TQwer (8 Panels) 


Joints of one of Dynamic analysis Static 

tbe leg members «— — ■ — * — — = — — — — analysis 

in eacb group 3 mode 6 mode 

superposition superposition 


2, 

8 

1 .518 

2.782 

7.125 

8, 

12 

4.850 

6.535 

10.166 

12, 

16 

14 .616 

15.560 

17.473 

16, 

22 

17.674 

I9i420 

22 .098 

22, 

26 

22.754 

25.415 

26,250 

26, 

30 

28.467 

34.235 

34.317 

30, 

36 

31 .775 

38 .504 

36.744 

36, 

40 

41 .087 

47.450 

45.267 

40, 

44 

44 .829 

49 .450 

48.370 

44, 

48 

49 .226 

50 .096 

52.013 

48, 

52 

54.790 

55 .111 

56,544 

(M 

56 

58 .642 

60 ,291 

58,817 

56, 

60 

62 .858 

65.848 

61 .594 

60, 

64* 

65.939 

*70,311 

63.159 

64, 

68 

69 .187 

74.782 

65.186 


Forces in metric ton 

* Extend tbe joint ntambering beyond 63 in the same sequence of 
Pig. 2.7. 
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Table 

5.6 Member Forces in 25m Body Tower (10 Panels') 






Joints of one of 
the leg members 
in each group 

Dynamic 

analysis 

Static 

analysis 

3 mode 

superposition 

6 mode 

superposition 

2, 

8 

1 .050 

2^536 

7.048 

8, 

12 

4.572 

6.552 

10,166 

12, 

16 

15.571 

17.174 

18.691 

16, 

22 

18.225 

21.044 

23.039 

22, 

26 

24,005 

27.088 

26.261 

26, 

30 

32 .214 

35.859 

34.308 

30, 

36 

33.953 

39 .458 

36.902 

36, 

40 

43.165 

50 .290 

45 .847 

40, 

44 

49.948 

56.142 

51.194 

44, 

48 

53‘.036 

57.623 

52 .428 

48, 

52 

59 .864 

62.336 

58 .200 

52, 

56 

65.138 

66.045 

61.115 

56, 

60 

71 .293 

72.118 

64.862 

60, 

64 

75 .819 

78-.651 

66 .842 

64, 

68 

80 .862 

85.669 

69.501 

68, 

72* 

84.607 

91.392 

71.011 

72, 

76 

88.552 

97.084 

73,051 


Fdroes in metric ton 

* Extend joint numbering similar to 20m body tower 
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5 *4 .2 Ea rthquake Load Analysis .Results 

Based on the recommendation £)f IS: 1893, wind load 
is not considered during earthquake. However other loads, 
viz,, dead loads due to the conductors and ground wire 
are considered. 

Ground acceleration due to earthquake is considered 
acting in the transverse direction (X direction) of the 
transmission line tower. Hence the ground acceleration in 
the other two directions, i,e., Y and Z directions, are 
zero. This is incorporated in the three dimensional 
dynamic analysis hy making necessary changes in the 
calculation of q. Eq, (5.41), In other words, the 

J 9 UlcUX. 

summation in Eq. (5.41) is carried out for i=1 ,4,7, . . ., (n£^2) 
only. 

In the eigenvalue analysis, we have observed that 
the lumped masses at joints are very small* This suggests 
that the loading due to earthquake ground motion would be 
small. Hence as a first attempt, analysis considering 
loads only due to ground motion is carried out. This study 
gives the stress in the members due to earthquake load 
only. Of course the stresses due to the other loads can 
be superimposed. 
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Maximum member forces due to earthquake load are 
computed by considering the contribution from the first 
6 modes. However, in this case the contribution due to 
the bending modes in the longitudinal direction as well 
as due to torsional modes is negligibly small. 

A tower of 25m body height with 10 panels in the 
body and base width of 4.2m is considered, Ihe member 
forces computed are given in Table 5 .7, assuming the 
tower to be located in the most critical earthquake zone. 

Table 5 .7 Member forces due to Earthquake Load 


Joint numbers* of one of 

the leg members in each Member force (ton) 

group (body portion) 


36, 40 
40, 44 
44, 48 
48, 52 
52, 56 
56, 60 
60, 64 
64, 68 
68, 72 
72, 76 


0.9271 
1 .0046 
0.9936 
1 .0298 
1 .0110 
1 .1050 
1 .2059 
1 .5100 
1 .3975 
1 .4824 


* Same as in Table 5 .6 
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Prom Table 5 .7, it is clear that the members are 
subjected to stresses of very low magnitude . Hence it is 
concluded that earthquake loads are not critical for the 
design of transmission line towers. 

5 *5 Re sults of Optimization and Discussion 

The optimum configuration design obtained under 
static loads and mentioned in Chapter 4 is taken as 
starting design for the present study of optimization 
under dynamic loads. Available I.S. angles are used for 
the members of the tower. Member selection is made to 
satisfy the design criteria as per I.S. 802-1973, discussed 
earlier in Section 4.5. The preassigned parameters for 
the optimum design under static load condition (page 82) 
are kept same for this study too . 

Optimum design is carried out for 15m, 20m, and 
25m heights of tower body. The results of optimum design 
are presented in Table 5 .8 and Table 5 .9 respectively for 
undamped system and for 2 percent damping. The CPU time 
taken for one function evaluation is also given in Table 5.8. 
The weight of the tower at the end of each iteration is 
plotted in Riga. 5 ,9 and 5.10 for undamped case and for 2 
percent damping respectively. 

From Table 5.8, it is observed that the optimum 
base width varies from 4,34m to 5.49m when the tower is 
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assumed undamped. However, the range of optimum hase 
width is 4.02m to 5.20m when 2 percent damping is considered, 
Table 5.9. Since 2 percent damping is recommended for 
steel structures, the latter values of optimum hase width 
are more realistic for transmission line towers , Contrary 
to the results of optimum design under static loads, there 
is considerable variation in the base width corresponding 
to the optimum tower imder d 3 mamic loads . Optimum base 
width increases for increasing height under dynamic loads 
while it remained more or less constant under static loads. 
Thus, for a more realistic optimum configuration design of 
transmission line towers, considerations in the dynamic 
response regime seem to be inevitable. 

Conclusions already made, from the study leading 
to optimum configuration under static loads, with respect 
to panel heights of transmission line towers remain unchanged 
under dynamic loading; i,e,, the number of panels should 
be chosen such that the height of central panelCs) is 
2 ,5m; the top and bottom most panel heights may be taken 
as 2,00m and 3.00m respectively; and the other panel heights 
may be fixed from linear interpolation of the top and bottom 
panel heights. 

The reduction in weight is maiximum corresponding 
to change in base width. This can be observed from Higs, 

5 .9 and 5 .10 wherein the reduction in weight in the first 



Weight, kg 


I 



Weight , kg 



iteration (wMch. corresponds to change in base width) is 
significant. The reduction in weight due to change in 
panel heights (other design variables in the present 
foimulation) is relatively small. 

The variation of in the process of optimization 
is plotted in Pig. 5. 11 (a) to (f) for the six eigenvalues 
considered in the present work. In each plot, the 
corresponding is plotted for undamped and with 2 percent 
damping of three towers considered. Prom these graphs it 
is observed that the change (increase) in the value of 
natural frequencies is maximum in the first iteration, 
i.e., corresponding to change in the base width of the 
tower, Purthermore, changes in are negligible for the 
corresponding changes in panel heights. In other words, 
small changes in panel heights resulting during process 
of optimization do not significantly affect the natural 
frequency. Since the change in panel heights during the 
process of optimization is small and the corresponding 
changes in natural frequencies of the tower are also 
negligible, the total computational work involved can be 
substantially reduced as follows; The optimum design 
considering the base width and panel heights may be carried 
out first for static loads and then one- dimensional 
minimization only with respect to the base width be 
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Fig5-11 (b) Variation of C 02 during optimization 
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undertaken in the dynamic response regime considering six 
eigenvalues from the lowest end of the spectium* 

In the present study, the change in natural 
frequency is an increase from the natural frequency of the 
corresponding starting design for all the transmission 
line towers considered, i.e., the tower becomes stiffen, 

This states that there is a reduction in the deflection 
at the joints. Hence it is justified from this observation 
that the constraints on deflection are not necessary for the 
optimum design of transmission line towers • 

The area of leg members in the starting design 
as well as in the optimum design is plotted in T'ig,5 .12(a) 
to (c) against midheight of panels above g 2 X)und. The 
reduction in the area of leg members is maximum in the 
body portion of the 25m body tower. In other words a 
considerable reduction takes place in the mass of the body 
of 25 m body tower. On the other hand, there is maximum 
increase in the base width of this tower, Table 5 .9 • Those 
two observations state that the optimum tower resist the 
dynamic loading by increasing the stiffness contribution 
and reducing the contribution due to inertia. 

A parametric study has been carried out to observe 
the effect of damping on the optimum design. For this 
purpose, range of damping of 0 to 5 percent is considered. 




starting design 

Optimum 
Damping-- 2% 

Height of body- 20m 


Area, cm 


Fig. 5*12 (b) Area of leg members along the height 
of tower 



Height , 


IG) 



Fig. 5-12 (c) Area of leg members along the height of tower 



The optimum weights obtained with different damping ratios 
are presented in Table 5*10. The corresponding optimum 
base widths are also given in this table. 

Prom the results of parametric study presented in 
Table 5.10, it is observed that there is a general tendency 
of reduction in weight with increased damping. However, 
the reduction in weight is not much. This is due to the 
fact that the difference in (Pig. 5.8) is almost 

negligible for the range of natural frequencies 
corresponding to the optimum design of transmission line 
towers. Prom the results presented in Table 5.10, it is 
further noted that the reduction in weight is not monotonic 
with increased damping. This is explainable again from 
the (I>LP)^g_y curves plotted in Pig. 5.8 which itself is 
not a strictly monotonic function of w for the ranges 
of w under discussion. Purthermore, the discrete choice 
of final available sections as members of the tower also 
affects monotonicity of the convergence. 
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COITCLUSIOITS MD SUGC-SSTIONS POR RmiHER WORK 

^ Introduction 

The configuration optimization of transmission line 
towers has been carried out both for static and dynamic 
behaviour of the tower. The configuration design variables 
considered are the base width of the tower and the pamel 
heights of tower body. The available angle sections have 
been used for the members of the tower. The resulting 
optimum configuration design problem turns out to be an 
unconstrained minimization problem which has been solved 
by using Powell *s method. 

The conclusions drawn from the results of the 
present study are grouped in Section 6,2, A set of 
recommendations for extension of the present work is 
given in Section 6.3, 

6 ,2 Conclusions 

(1) It is sufficient to model the transmission line 
towers as a space truss both for static and dynamic analysis. 
This conclusion for dynamic loads is based on the comparison 
of the eigenvalues of the towers and that of the individual 
members of the corresponding tower. This study establishes 
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that the response of the overall tower structure precedes 
that of individual members in the present case because the 
individual members are constrained to be designed for a 
l/r ratio of maximum of 200. This maintains the natural 
frequency of longest individual member by and large constant 
for towers of all sizes where as the frequency of the tower 
goes on decreasing for increasing tower height, 

(2) The three critical loading conditions on the tower 
are the broken conditions of the ground wire, top conductor, 
and the middle conductor which need be considered for the 
purposes of design of transmission line tov/ers. Earthquake 
loads are not critical for the design of transmission line 
towers . 

(3) Stiffness approach is the most efficient method 
for the analysis of transmission line towers from the 
computational view point . 

(4) Inorder to have a realistic response prediction 
for wind load, considerations of equivalent static loads 

are insufficient and a rigorous dynamic analysis is inevitable. 
This is more so when optimum design is the goal, 

(5) A check, particularly for small towers, for safe 
response is warranted by the study of eigenvalues of the 
overall tower and of the longest member therein. This is 
due to the fact that in normal practice a conservative design 
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for small towers shall he produced by providing a wide 
base resulting in large stiffness of the tower, Ihis 
may result in the overall natural frequencies of the 
tower higher than the natural frequency of the longest 
member in such towers. In that situation the dynamic 
response of the longest member shall govern the safe 
design . 

(6) Transmission line towers are structurally quite 
stable and hence it is not necessary to embed overall 
stability considerations in the optimum design process, 

(7) The variation in base width has marked effect 
on the weight of the tower. The optimum weight is more 
sensitive to this design variable in the dynamic response 
regime as compared to the static behaviour, Prom the 
results, it is observed that the base widths corresponding 
to optimum design tower for 15m, 20m, and 25m body heights 
are around 4.00m, 4.50m, and 5.40m respectively. In the 
optimum design under static loads these values turned out 
to be around 4,00m for all towers considered. The latter 
phenomenon in itself is not explainable. However, it gets 
corrected in the dynamic response regime, 

(8) The number of panels in the body of the tower 
should be chosen such that the height of central panel (s) 
of the body is around 2,5m, Further, the top most panel 
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height may be kept around 2,00m and that of bottom most 
panel around 3.00m for achieving minimum weight tower. 

The intermediate panel heights may be fixed by linear 
interpolation of heights of top and bottom panel. 

( 9 ) The panel heights corresponding to optimum design 
of the tower turns out to be the same both under static and 
dynamic behaviour. Hence the panel heights may be decided 
from the optimum design under static loads , This leaves 
only a one- dimensional optimization to be performed with 
respect to the base width under dynamic loads so that 
overall computational time is minimized and realistic 
optimimi design of transmission tower is achieved. 

(10) The reduction in weight of tower for increase in 
damping is negligible. 

(11) The conventional model of plane frame for the 
dynamic analysis of slender towers is confirmed to be 
reasonably accurate through the present study. However, 
since slenderness is only a relative concept, if the loading 
is such as to cause rotation about the vertical axis, then 
spflCe truss model for dynamic analysis as concluded in 

( 1 ) above will have to be adopted , 

(12) The available Indian Standard angle sections are not 
rational, particularly the heavier ones from the view point 
of optimum design. These can be revised (proposal presented 
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in Appendix 2) so as to make them more rational from the 
point of view of minimum weight design of such structures. 

6 SQjop.e. f03? jPu rther lork 

The present work has been oarried out for 220 K7 
transmission line towers. The same can be extended for 
440 KY and higher transmission voltages . This may require 
towers of very different configurations , 

Wind loading is a random phenomenon. Therefore, it 
should be treated as random load to get more accurate 
response calculations and hence more rational optimum 
design. The random vibration analysis, of course , requires 
statistical data on wind over a reasonable period of time 
which is to be collected. 

It has been observed that broken wire conditions 
occur very rarely, and in design, these are the critical 
loads. Hence a study of the factors of safety is essential, 
furthermore, it would be more economical to design the 
tower without consideration of the brokenwire conditions 
and at the same time to prevent the total collapse of the 
tower in case a conductor breaks. A possible solution to 
this is to design the cross arms such that the cross arm 
turns away from the tower under the unbalanced tension in 
the event of breaking of a conductor. The rest of the 
tower structure would then be safe. Even material having 
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APPENDIX 2 


RATI ON AC MG-LE SECTIONS 

I ntroduction. 

jingle sections are extensively used as tension or 
compression members in industrial "buildings, transmission 
towers, trusses etc. These widely used angles should have 
optimum dimensions so that their use as truss member will he 
economical^^*^ , 

Most of the angle sections that are listed in 
ISI Handhooh^^ are found to he uneconomical. This observation 
is brought out in this study and' the proposed rational sections 
are listed. Rational angles are derived with the aim that 
these angles are used mainly for truss members . The study 
is restricted to equal angle sections only. 

S trength of Column Members 

The maximum load that a column can carry without 

buckling is given by^^ 
ti^EI 

P * c (i2.1) 

1 ^ 

where Q = constant which depends on the boundsiry conditions 
at the two ends (=1 for hinged ends, 4 for built- 
in ends) , 
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E = modulus of Elasticity of the material, 

I = least moment of inertia, i.e,, minimum of the 
two principal moments of inertia, 
and 1 = effective length of columri. 

From the above expression, it can be seen that the load 
carrying capacity of a column of given effective length is 
directly proportional to the least moment of inertia. That 
is, the value of I has to be maximized by properly dimensioning 
the cross section, to get maximum permissible compressive load. 
The relationship between moment of inertia and area of the 
section, 1, is given by 

I = A r^ (i2.2) 

where r = radius of gyration. 

To achieve maximum economy, the section chosen 
should have minimum area and maximum radius of gyration. 

In other words maximize the least moment of inertia of the 
section for a given cross sectional area'. This is possible 
by suitably dimensioning the section. 

Sectional Properties of the Available Angles 

To study the properties of the available angles, the 
least moment of inertia of these angles are listed against 
their increasing area of cross section is- Table A2%1 . These 
values are also plotted in Fig,A2.1. From Table A2.1* it can 



Table A2,1 Properties of the Indian Standard Angle Sections 


Moment of Inertia for Increasing 4rea of Angles 


SI; 

Mo. 

Area 

cm^ 

^min 

cm 4 

SI. 

Mo. 

Area 

cm 2 

I . 
min 

cm4 ^ 

SI. 

Mo. 

Area 

cm 2 

%in 

cm 4 

1 

1*12 

0,2* 

25 

5.68 

5.3 

49 

14.02 

29.4 

2 

1 .41 

0.3* : 

26 

5.75 

7.7^ 

50 

15.05 

36.0 

3 

1.45 

0.2 

27 

6.25 

9.9* 

51 

15.-39 

58.4* 

4 

1.73 

0.6* 

28 

6.26 

7.0 ^ 

52 

i7;o2 

78.2* 

5 

1,84 

0.4 

29 

6.77 

12.5* 

53 

17.03 

51.6 

6 

2,03 

0.9* 

30 

6.84 

9.1 ^ 

54 

17.81 

42.4 

7 

2.25 

0.5 

31 

7.27 

15.5* 

I 55 

19.03 

71*8 

8 

2.26 

0.7 

32 

7.44 

11.7 

! 56 

20.19 

60.9 

9 

2.34 

1.4^ 

33 

8.06 

14.8 

I 57 

20.22 

131.4* 

10 

2.64 

2.0* 

34 

8,18 

9.1 

^ 58 

21.06 

96.3 

11 

2.66 

1.2 

35 

8,66 

18.4* 

59 

22.59 

84.7 

12 

2.77 

0.9 

36 

8.96 

11.9 

60 

25.02 

113.8 

13 

2.95 

2.8* 

; 37 

9.29 

22,5* 

61 

25.06 

162.1* 

U 

3.07 

1.8 

38 

9.76 

15.3 

62 

29.03 

249.4* 

15 

3.27 

1.5 

39 

10.02 

11.2 

63 

29.82 

191.8 

16 

3.47 

2.6 

40 

10.47 

32 , 0 * 

64 

30.81 

139.5 

17 

3.78 

2.2 

41 

10.58 

19.3 

65 

54.59 

296,0# 

18 

3.86 

1.7 

42 

11.00 

14.6 

66 

36.81 

. 235,0 

19 

3.88 

3.6* 

. 43 

11.38 

24.0 

67 

42,78 

36-3,8* 

20 

4.28 

3.2 

44 

11.67 

44.5* 

68 

46.61 

715,9* 

21 

4.47 

2.6 

45 

12.00 

18.8 

69 

50.79 

429.5 

22 

4.79 

4.5* 

46 

12.21 

29.4 

70 

57.80 

883.7* 

23 

5.07 

3.8 

47 

13.02 

23.7 

71 

68,81 

1046.5* 

24 

5.27 

5.9* 

48 

13.79 

42.0 

72 

93.80 

1411*6* 
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be observed that some of the sections (marked #)have special 
sectional properties, among the available sections , The 
least moment of ineirfcia of these sections is increasing 
monotonically as the area increases. There are 29 sections 
out of the total 72, having such special property. Hence 
these sections will be preferred to others in the design 
of struts. For example , the section having area 2.95 cm^ 
would be recommended instead of 3*86 cm^ since the least 
moment of inertia of the latter is less than that of the 
former, i.e., load carrying capacity of the former is greater. 

From Fig*A2»1, it can be observed that the above 
special sections fall far away on the left from the mean curve 
The equation of the mean curve is obtained by least square 
technique and the relationship between A and l is given by, 

A= 2.48294 (l)°«^'7276 

In an automated design, using the sections listed 
in Table A2.1, the sections picked up would be one of the 
angles marked * in Table A2,1, if the only consideration in 
the design is Euler buckling load. Other soctionsj if 
adopted will lead to uneconomical design. 

Permissible Stresses in .Angle Sections 

From the dimensions of the 29 special sections 
(marked * in Table A2.1), it can be observed that the b/t 



mm , 



.Ftg.A2.1 Variation of moment of inertia with area 
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ratio (defined below) varies from 13 to 16 for most of these 
sections. In general the least moment of inertia can be 
increased infinitely by increasing the b/t ratio for a given 
value of area. But the allowable stress in such a section 
gets reduced so as to avoid local buclcLing of such thin 
flanges . 


IS:802 (part l)'-1973^^ stipulates reduction in 
permissible stress as the b/t ratio increases. For steel 
confirming to IS ;226>-l969® these formulae are given as 
follows : 

( i) For the largest effective slenderness ratio(— ) < 120; 

(Su)^ 

^ ? 

) 2600 kgf/cm where b/t <13 

(~; V o D 


( 1.0 - 


Ji 


F^ = (1.0 - 

a 31200 

Kl 2 
(— ) ^ 

^a ■ ” 3^200 ■ 

KI 

(ii) For ( ) > 120 

r 


■) (4680-160 •—) kgf/cm 13 < 


t “ 


590000 

(b/t)‘ 


( A2 .4t>) 


for b/t > 20 
(A2.4c) 


20,000,000 , 2 

j = ^ J kgf/cm 

(S)2 

r 


(42.5) 


where Fg^ = permissible stress in compression 
K = fixity factor 

L = length of the compression member 


20 


@ Specification for structural steel (standard quality) 
fourth revision. 
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KL 

= largest effective slenderness ratio of any 

unbraced segment of the member 

b = distance from the edge of fillet to the extreme 
fibre 

and t = thicloiess of material, 

KL 

When > 120, the pemissible stress, Eq.(A2.5) 

is independent of (b/t) ratio and it is nothing but Euler 
buckling load, Eq., (A2,1) divided by area, i,e., when the 
member is slender (b/t) ratio can be large since the permissible 
load is Euler buckling load itself. In practice, slender 
members are adopted only for secondary bracings and for few- 
primary diagonals. Hence, for these cases, any one of the 
nominal sections can be used. 

Most of the leg members of the transmission towers 
and main members of the trusses have to be designed for heavy 

axial loads and also they are not slender. Since these members 

\ 

of the structure form the major portion of the materials 

consumed, they should be designed economically to gain overall 

reduction in the materials. Hence attention has to be 

,KI 

focussed on the members having (~) ratio less than 120, 

Rational ingles 

With the aim of reducing material consumption, one 
would, naturally like to take the angles having minimum area. 



1,88 


^ea of the section, will he least when the allowable stress 
of the section is maximum. Purthemore, as discussed in 
the earlier sei^tions, the dimensions of the angle should be 
chosen with the aim of having minimum area and maximum strength. 


let Bq be the distance from the edge of one flange of 
equal angle (Big. A2.1) to the centre line of the other and 
k is the b/t ratio. 



Area of cross section, 



A=2Bot 


( 2k4- 1 ) p 
(2k+1) 


and 


dA 8Bp 

die ” ~ (2k+1) 


(A2.6a) 

(A2.6b) 


Brom the ISI Handboolc?^, it is seen that the radius 

of gyration is almost constant for different thickness, t, of 

angles of a given flange width B and hence it is independent of 

thickness t. This can also be observed from Big,A2.1, where 

the Imin versus A is almost a straight line for a fixed flange 

width. So the quantity I 1 ,0 - ( ^ Bqs,(A2.4) 

- 31200 r 

is a constant for a fixed value of B. 



(A2.7) 
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and 


dp dA dPa 

dk ~ dk " a ^ 


Substituting the expressions from Bq, (12.4) 
dp dA 

dk dk 

dA 
dk 


dA 


2600 C, k < 13 (12 .8a) 

[^2600-160 (k- 13 )Jo V. Ax: 160 C;13 < k < 

r dA -1 

2600 G - — (k-13) + A / 160 G 

t dk 


dk 

dA 590000 


where G =| 1 .0 - 3^200 


dk k^ 

L_ ( 7 

rsr^r^ -v* ^ J 


590000 x: 2 
C- A 0 


k5 


( 12 .8h) 

\ 

20 

(A2.8C) 


The above equations reveal that for unit change in area of 
the section the corresponding change in capacity of section 

is maximum in the sections for which k _< 13. Hence sections 

’ b -b 

having ~ ^ 15 are well utilized than those having >13. 

In otherwords, the sections for which the allowable stress 

is maximum are rational. In the earlier sections it was 
b 

observed that ratio can be chosen as high as possible to 

get better sectional properties. These observations point 

to the fact that the sections should be dimensioned in such 
b 

a way that ^ 13 to get optimum use. 

Recommended Angle Sections 


20 


Based on the observations in the previous sections 
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the proposed rational angles are listed in Table i2.2. Is 
far as possible the proposed sections have k value equal to 13. 
The sectional properties of these angles are plotted in 
Pig. A2 .1 along with the mean curve for the presently 
available angles . 

When the sections are chosen to have rational dimensions, 
the number of sections obtained are a few for the standard 
width of angles available in ISI Handbook, So additional 
sections having width in between those available are also 
recommended, while doing so, only those angle sections having 
monotonic increase in sectional properties are presented. 

There are 48 such sections, listed in Table a 2 ,2 , These are 
closely spaced. Hence these angles when used in an automated 
design would produce rational design in comparison to the 
presently available sections . 

The equation of the mean curve for the recommended 
angles in terms of A and is given by 

A= (A2.9) 

The deviation of the sectional properties of the proposed 
angles from the above mean curve (not drawn in Pig.A2.1) is 
very small. Hence, in case a continuous function in terms 
of area and moment of inertia is required, in any optimization 
study, Bq, (A2.9) could be used without much of error. 
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Table .12, 2 Proposed Hqual An.g:le Sections 


Sl.l'Io, 

Siae 

AxiBxt 

mm 

Area 

cm2 

^min 

cm4 

1 

20x20x3 

1.12 

0.2 

2 

20x20x4 

1.45 

0.2 

5 

25x25x3 

1.41 

0.3 

4 

25x25x4 

1.84 

0.4 

5 

30x30x3 

1.73 

0.6 

6 

30x30x4 

2.26 

0.7 

7 

35x35x3 

2.03 

0.9 

8 

35x35x4 

2i66 

1.2 

9 

40x40x3 

2.34 

1.4 

10 

40x40x4 

3.07 

1 .8 

11 

45x45x3 

2.64 

2.0 

12 

45x45x4 

3.47 

2.6 

15 

50x50x4 

3.88 

3.6 

14 

55x55x4 

4. 26 

4.9 

15 

55x55x5 

5.2? 

5.9 

16 

60x60x5 

5.75 

7.7 

17 

65x65x5 

6.25 

9.9 

18 

70x70x5 

6.77 

12.5 

19 

75x75x5 

7.27 

15.5 

20 

75x75x6 

8,66 

18.4 

21 

80x80x6 

9.29 

22.5 

22 

85x85x6 

9.87 

27 i6 

23 

90x90x6 

10.47 

32.0 

24 

90x90x7 

12.14 

37.2 

25 

95x95x7 

12.84 

44.8 

26 

100x100x7 

13i54 

52*6 
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Table A2,2 contd,... 


31. No. 

Size 

AxBxt 

mm 

Area 

2 

cm^ 

^min 

cm4 

27 

1 05x1 05x7 

14.25 

61.3 

28 

1 05x1 05x8 

16.20 

69.2 

29 

1 1 0x1 1 0x8 

17.02 

78.2 

30 

1 20x1 20x8 

18.62 

104.2 

31 

120x120x9 

20.85 

116.6 

32 

1 30x1 30x9 

22.65 

146.0 

33 

130x130x10 

25.06 

162.1 

34 

1 40x1 40x1 0 

27.05 

206.2 

35 

150x150x10 

29.03 

249.4 

36 

150x150x11 

3t.82 

277.6 

37 

1 60x160x1 1 

34.08 

340.0 

38 

160x160x12 

37.05 

368.2 

39 

170x170x12 

39.45 

' 444.4 

40 

1 80x1 80x1 2 

41.83 

528.8 

41 

180x180x13 

45.18 

569.0 

42 

190x190x13 

47.78 

672.9 

43 

190x190x14 

51.31 

720.1 

44 

200 x 200 x 1 4 

54.09 

841.9 

45 

200x200x1 5 

57.80 

883.7 

46 

200x200x16 

61.49 

95U2 

47 

200x200x17 

65^.16 

1005,0 

48 

200x200x18 

68.81 

1046.5 
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The use of the proposed angles in transmission line 
tower design has been described in Section 4.5. The optimum 
weights of tower obtained with the proposed angles are 
compared to the one obtained by using the available sections, 
Table 4.6. I’rom the results of the optimization study, it is 
confirmed that the recommended sections (Table a 2 .2) could 
be used economically in the construction of any truss-like 
structures . 



APPENDIX 3 


RAYLEIGH- RIIZ MEIHOD 


The subspaoe iteration method uses the Rayleigh- Ritz 
principles and the Eq« (5.10) to Eq. (5.14) are nothing hut 
the steps involved in Rayleigh- Ritz method. Hence the 
Rayleigh- Ritz method is briefly described below: 


let denote the n-dimensional space in which the 
operators K and M are defined. The Rayleigh minimum principle 
states that 


X ^ = min. P (^) 


(A3.1) 


in which the minimtun is taken over all vectors j6 in Y^ and 
P(0) 

is the Rayleigh quotient defined as 


P(B) = 


->T 

^ K ^ 

M ^ 


(A3 .2) 


The nmerator in the above equation is twice the maximum 

strain energy for a given ^ (displacement vector in 

structural application) and the denominator is twice the 

2 

maximum kinetic energy of the structure, divided by w , 

The minimum of p(0) can be obtained by any one of the methods 

Co 

of unconstrained minimization « 

In the Ritz analysis all the vectors p in a q- 
dimensional subspace Vg^ of Y^ are considered, A typical 





element in the subspace is given by linear combinations of 
Ritz basis vectors > 




i=1 


^1 ^1 


(A5.3) 


in which a^- = the Ritz coordinates. Substituting Sq, (A5.5) 


in Eq« (A3. 2), q_ 


P(j3) = i 


2 S a . a. k. . 
=li=1 D 1 ID 


1 q 

E 2 a a. m . . 
D=1 i=1 ^ 


(A3.4) 


-T - 

is obtained with k. ^ = Xj^ K x 


and 


m- . = X. 
ID 1 


(A3.5) 




3p(^) 


The necessary condition for a minimum of P (p) is g-™ = 0, 

i=1,2, ,,.,q, because a^ are the only variables. This gives 
K a =~^ M a (A3. 6) 


where a is a vector listing the Ritz coordinates; K and M are 
matrices of size q x q with typical elements defined in 
Bq. (A3. 5)* 

The solution to Eq, (A3, 6) yields q values“V ^ 2» • * • 

and the corresponding vectors , ....... pq > which are 

obtained using Eq. (A3. 3)* The values X^, i=1,....,q are 
97 98 

upper bound ^ approximations to the exact eigenvalues of 
the original problem, Eq. (5*4)> i.e., 



(A3. 7) 


1 - 1 


^2 — ”^ 2 ’ 


.2 

q. - 


W- <1 


q. 


Tile error in the solution depends on the Ritz basis vectors 
chosen because the approximate eigenvectors , . . . . , 
the elements of the subspace. 


"^1 1 



iPPEJ^DIX 4 


SOLUTlOi4 Of ;ERALIZED SYMI^TRIQ EIG-ENViilillS PROBLEM 

'LLe subspace iteration method requires the solution of 

the generalized symmetric eigenvalue problei^ Eq, (5.13). 

Since the size of the matrices involved are only qxq and also 

it is necessary to calculate all the q eigenpairs, Jacobi 

84 

method recommends itself . 

The method of Jacobi to solve the generalized eigen- 

problem 

(A - XB ) g = 0 (A4.1) 

is presented, here, with matrices A and B symmetric and 
known. In the presentation below, A will be same as and 

B is same as in Eq. (5.13). Phe requirement in this 

method is that B should be positive definite. In the present 

-T 

case, B = M^^^= M is nonsingular, since it is 

quadratic* 

Eirstj the positive definite matrix B of the generalized 
symmetric eigenvalue problem (A - XB) g = 0 is to be transformed 
by Jacobi's method to a diagonal form D, using the orthogonal 
matrix E: 


T 

D = E B E 


(A4.2) 



1.98 


Tlie generalized eigenproblem is transformed with the help of 

matrix E in that, on one hand, the identity matrix is 

CD 

introduced in the form EE , and on the other, the eigenvalue 

T 

equation is multiplied from the left by E : 

E^ (A - XB) EE^ g = 0 (A4.3) 

Eq, (A 4 . 3 ) is rewritten, using Eq. (A4,2) as 

(E^ A E - XD ) E^ g = 0 (A 4 . 4 ) 

CDhe generalized eigenvalue problem is hereby very nearly 
converted into the special one. In place of the identity 
matrix I, Eq. (A4,4} has a diagonal matrix D with only 
positive elements which ane equal to the eigenvalues of B, 

CDhe diagonal matrix consisting of the square roots of the 
positive diagonal elements of D, exists, as does its inverse 
Thus Eq. (A4.4) can be further transformed into 

E"^(E^ a E - XD) if (E^ g) = 0 
or (D'^’e^ a B XI)(B^ g) = 0 (A4.5) 

i m i 

The matrix T = B A E D"" (A4.*6) 

must be symmetric since A is S3mmietric, E is orthogonal, 

— ‘i’ 

and B is a diagonal matrix. Through 

h = B^ E^ g (A4.7) 

the generalized eigenvalue problem is reduced to a special 
one with symmetrix matrix T, i.e,, 
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(T - X I ) h = 0 


(A4.8) 


The eigenTalues X of T agree vath that of the 
generalized eigenvalue problem (A - XB) g = 0 , and the 
interrelationship of the eigenvectors of h and g are given 
by (A4-.7}. The special eigenvalue problem, Eq,(A4.8) can be 
again solved by Jacobi's method. Here we encounter the 
transformation matrix Y containing the eigenvectors of T 
in its columns. Through Eq, (A4.8) the matrix G- of the 
eigenvectors of the problem Eq, (A4,1) can be obtained. 


The algorithm for finding the eigenvalues and 
eigenvectors g^ of a generalized S3nmiietrix eigenvalue problem 
(A-^^ B) gj^ = 0 consists of the following stepsJ 


,T 


(a) E B E = D 


(Jacobi) 


-2 ^T 


(b) T = I)“ E A E D 




(c) 


(A4.9) 


T Y = D 


(d) G = E D‘ 




1 


(Jacobi) 


Eigenvalues X^^ are the diagonal elements of and the 

corresponding eigenvectors are the appropriate columns of 

matrix G, Matrix G is itself not orthogonal, instead we 
T 

have G B G = I, where the generalized metric is expressed, 


The method of Jacobi, that is required in (A4,9) 

84,90,96 


is well documented in literature 
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